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ABSTRACT 

We present the first high spatial resolution X-ray study of the massive star 
forming region NGC 6357, obtained in a 38 ks Chandra/ ACIS observation. In- 
side the brightest constituent of this large HII region complex is the massive 
open cluster Pismis 24. It contains two of the brightest and bluest stars known, 
yet remains poorly studied; only a handful of optically bright stellar members 
have been identified. We investigate the cluster extent and Initial Mass Function 
and detect ~800 X-ray sources with a limiting sensitivity of ~ 10^° ergs s~^; 
this provides the first reliable probe of the rich intermediate-mass and low-mass 
population of this massive cluster, increasing the number of known members 
from optical study by a factor of ~ 50. The high luminosity end (log [2-8 
keV]> 30.3 ergs s^^) of the observed X-ray luminosity function in NGC 6357 is 
clearly consistent with a power law relation as seen in the Orion Nebula Cluster 
and Cepheus B, yielding the first estimate of NGC 6357's total cluster popula- 
tion, a few times the known Orion population. We investigate the structure of 
the cluster, finding small-scale substructures superposed on a spherical cluster 
with 6 pc extent, and discuss its relationship to the nebular morphology. The 
long-standing Lx — 10~^ Lf,oi correlation for O stars is confirmed. Twenty-four 
candidate O stars and one possible new obscured massive YSO or Wolf-Rayet 
star are presented. Many cluster members are estimated to be intermediate- 
mass stars from available infrared photometry (assuming an age of ~ 1 Myr), 
but only a few exhibit i^-band excess. We report the first detection of X-ray 
emission from an Evaporating Gaseous Globule at the tip of a molecular pillar; 
this source is likely a B0-B2 protostar. 

Subject headings: ISM: individual (NGC 6357) - open clusters and associations: indi- 
vidual (Pismis 24) - stars: formation - stars: mass function - stars: pre-main sequence 
- X-Rays: stars 
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1. Introduction 

NGC 6357 (= W 22 = RCW 131 = 
S 11) is a large Hll region complex in 
the southern sky excited by the OB as- 
sociation Pismis 24, at a distance of ~ 2 
kpc (Neckel 1978; Felli et al. 1990; Massi 
et al. 1997; Bohigas et al. 2004). This 
young stellar cluster and its surrounding 
gaseous environment are relatively poorly 
studied; only a handful of stellar mem- 
bers have been identified and character- 
ized, and the molecular cloud has not been 
well-mapped. In optical and near-infrared 
images, the exciting cluster is not promi- 
nent against the Galactic field population 
except for a few bright members (Figure 
1). Yet, two of these stars are remark- 
ably luminous and massive with spectral 
types 03 If (M~ 200-300 Mq) and 03 III 
(M~ 100 Mo; Massey et al. 2001; Walborn 
et al. 2002). Such massive stars are usu- 
ally found in only the richest Galactic clus- 
ters and the giant starburst regions such 
as Carina, 30 Doradus, and NGC 3603. 
This raises the question is Pismis 24 much 
richer than it appears, or is it a poor clus- 
ter with a top-heavy Initial Mass Function 
(IMF)? This would be highly unusual, as 
top-heavy IMFs are seen only in extreme 
starburst environments such as the super- 
star clusters of M 82 (Elmegreen 2005). 

We investigate this question by measur- 
ing, for the first time, the low-mass pre- 
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main sequence (pre-MS) population of the 
Pismis 24 cluster and its environs. This 
is achieved with a sensitive X-ray obser- 
vation with the Chandra X-ray Observa- 
tory. Chandra observations can penetrate 
heavy absorption up to Ay — 500 (Grosso 
et al. 2005) and resolve point sources with 
sub-arcsecond resolution. In most cases, 
the X-ray detected stars appear on existing 
near-infrared (NIR) images; the advantage 
is that the strong X-ray emission of pre-MS 
stars (elevated 10^ — 10"^ above MS stars 
Preibisch & Feigelson 2005) traces mag- 
netic activity rather than photospheric or 
circumstellar disk emission, therefore effec- 
tively discriminating cluster members from 
older Galactic field stars. We estimate 
that the cluster has ~ 10, 000 stars, ~ 5 
times richer than the known Orion Neb- 
ula Cluster. The presence of two 03 stars 
is compatible with a cluster of this popu- 
lation. The known massive stars and the 
large number of low mass X-ray stars sug- 
gest that Pismis 24 is a very rich cluster 
with a standard IMF. 

Our study also gives considerable new 
information about the cluster and its X- 

ray properties. We provide a census of 
~750 cluster members with sub-arcsecond 
positions, 613 of which have NIR counter- 
parts. Although not complete, the survey 
includes pre-MS stars with masses extend- 
ing down to ~ O.3M0; most do not have 
circumstellar disks revealed by A'-band ex- 
cess. We study the spatial distribution of 
the cluster members and identify a possi- 
ble subcluster. We examine the ringhke 
morphology of the NGC 6357 nebula and 
investigate the possible feedback from the 
massive stars to their environs. We mea- 
sure X-ray properties of a dozen known O 
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stars, and list over 20 new candidate O 
stars. A faint X-ray source is found asso- 
ciated with a protostar in an Evaporating 
Gaseous Globule (EGG) for the first time. 
A few pre-MS stars are seen with very pow- 
erful X-ray flares with peak ~ 10^^ ergs 

Section 2 reviews past study of the Pis- 
mis 24 chister and NGC 6357 region. The 
Chandra observation and its analysis are 
described in §3, and a NIR observation us- 
ing the Simultaneous 3 color InfraRed Im- 
ager for Unbiased Surveys (SIRIUS) cam- 
era on the Infrared Survey Facility (IRSF) 
telescope is summarized in §3.5. The as- 
sociation between X-ray sources and op- 
tical and near-infrared (ONIR) stars, as 
well as the infrared colors of X-ray coun- 
terparts are presented in §4. The cluster 
population inferred from the X-ray emit- 
ting stars, their X-ray luminosity function 
(XLF), the spatial distribution of the clus- 
ter members, and the morphology of the 
nebula are discussed in §5. X-ray proper- 
ties of O stars and intermediate- to low- 
mass pre-MS stars are presented in §6. 
The possibility of diffuse X-ray emission 
is examined in §7. A brief summary of 
our findings is given in §8. This study 
is part of a series of X-ray investigations 
of high-mass star formation regions that 
includes the Rosette Nebula (Townsley et 
al. 2003; Wang et al., in preparation), M 
17 (Townsley et al. 2003; Broos et al., in 
preparation), the Orion Nebula (Feigelson 
et al. 2003; Getman et al. 2005 and asso- 
ciated articles), 30 Doradus (Townsley et 
al. 2006a, 2006b), and Cepheus B/0B3b 
(Getman et al. 2006). 



2. Observational Overview of NGC 
6357 and Pismis 24 

NGC 6357, a large H ll region complex 
showing an annular morphology in the ra- 
dio and optical (Haynes et al. 1978; Lortet 
et al. 1984), is located in the Sagittar- 
ius spiral arm and spans 60' x 40' in the 
southern sky. NGC 6334, another promi- 
nent star- forming complex likely associ- 
ated with NGC 6357, lies ~ 60' away 
(Neckel 1978; Ezoe et al. 2006). The mas- 
sive open cluster Pismis 24 (= C 1722- 
343 = OCl 1016) lies in NGC 6357's cen- 
tral cavity (Pismis 1959; Alter et al. 1970). 
The kinematic distance to NGC 6357, d = 
1.0 ± 2.3 kpc, is given by Wilson et al. 
(1970). Neckel (1978) derives d = 1.74 ± 
0.31 kpc with improved accuracy based on 
UBV and H/3 data. Spectroscopic paral- 
lax of 10 high-mass members in Pismis 24 
gives a larger distance of 2.56 ± 0.10 kpc 
(Massey et al. 2001), which we adopt in 
this paper. We should note that recent 
work (Kharchenko et al. 2005; Prisinzano 
et al. 2005; Arias et al. 2006) has found 
smaller distances to other high mass star 
forming regions (e.g., M8) than previous 
studies; we will discuss the uncertainties 
in distance-related stellar properties wher- 
ever applicable. 

Figure lb shows the large-scale mor- 
phology in the mid-infrared, observed with 
the Midcourse Space Experiment (MSX) 
satellite. A large ring-like structure ex- 
tends ~ 60' in diameter and bright nebu- 
losity is seen at the northern rim. The cen- 
tral cavity appears free of large amounts 
of warm dust. Far-infrared (FIR) contin- 
uum revealed several luminous (L~ 10^ 
L0) embedded sources coincident with CO 
line and radio continuum emission peaks 
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(McBreen et al. 1983). G353.2+0.9, which 

we study here, is the youngest {t ~ 10^ 
yr) and brightest of NGC 6357's three Hll 
regions. It hes at the northern interface be- 
tween the cluster and molecular cloud and 
contains a number of embedded sources 
(Figure la; Frogel & Persson 1974; Sakel- 
lis et al. 1984; Persi et al. 1986; Massi ct 
al. 1997; Bohigas et al. 2004). In contrast 
with NGC 6334, almost no water maser 
emission is found in NGC 6357 (Healy et 
al. 2004a). This may either indicate that 
massive star formation has ceased in the re- 
gion (Persi et al. 1986) or that it has been 
overrun by the ionization front (Healy et 
al. 2004a). G353. 1+0.6 is a more evolved 
Hll region viewed cdgc-on, bounded on 
the northern side by a different molecu- 
lar cloud. It contains a single 05 star, 
several B stars and embedded IR sources, 
but no indication of current star forma- 
tion can be found (Felli et al. 1990; Massi 
et al. 1997). G353.2+0.7, which is also 
a FIR peak, does not contain early type 
stars (FeUi et al. 1990; Massi et al. 1997), 
but it hosts the only water maser found 
in NGC 6357 (Sakelhs ct al. 1984). How- 
ever, Persi et al. (1986) identifies this water 
maser source as an evolved late M star. 

Despite a number of optical, infrared, 
and radio studies on NGC 6357, no X- 
ray observation of this field has been re- 
ported. Figure If shows an unpublished 9 
ks ROSAT PSPC observation revealing a 
few strong X-ray sources, including NGC 
6357's three H ll regions. The brightest X- 
ray emission coincides with the core of the 
Pismis 24 cluster, but the 05*74 T instru- 
mentation was unable to resolve the source 
population. 

The OB cluster Pismis 24 is centered ~ 



1' south of the G353.2+0.9 ionization front. 
Early optical study of NGC 6357 and Pis- 
mis 24 revealed ~ 20 0-type and early 
B-type stars (Moffat & Vogt 1973; Neckel 
1978, 1984; Lortet et al. 1984), including 
the unusally late-type WC7+07-9 Wolf- 
Rayet (WR)/0 star binary HD 157504 
(= WR 93; van der Hucht 2001). Two 
of the cluster members, namely Pismis 
24-1 (=HDE 319718) and Pismis 24-17, 
were recently classified as spectral type 
03.5, some of the brightest and bluest stars 
known (Massey et al. 2001; Walborn et al. 
2002). Only a dozen other 03 stars are 
known in the Galaxy (Massey & Johnson 
1993; Walborn 1994; Mafz-Apellaniz et al. 
2004). The two luminous 03 stars are eas- 
ily identified close to the sharp boundary 
of the bright optical nebula in the opti- 
cal image, accompanied by several other 
bright cluster members extending ~ 4' 
north-south (Figure le), but are not as ap- 
parent in the Two Micron All Sky Sur- 
vey (2MASS) K-band image (Figure Ic). 
Figure Id shows a SIRIUS JHK com- 
posite image, which has higher resolution 
and penetrates deeper into the cloud than 
2MASS. The two luminous O stars are ev- 
ident in the center of the SIRIUS image, 
which contains thousands of stars. Other 
noticeable features are the bright nebula 
with tenuous diffuse filaments in the north 
and the infrared dark column occupying 
the south-east corner of the image, pre- 
viously noted as the south eastern com- 
plex (S.E.C.) in CO mapping (Massi et 
al. 1997). Several heavily reddened stars 
can be seen in the S.E.C. region. However, 
most NIR stars are foreground and back- 
ground objects, not cluster members. 

Although most researchers believe G353.2-I-0.9 
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is excited by the hot stars in Pismis 24 
(e.g., Lortet ct al. 1984), FelU et al. (1990) 
argued that Pismis 24 is projected onto 
the H II region by chance, and the embed- 
ded massive stars (IRS 1, 2, and 4; Persi 
et aL 1986) supply the ionizing radiation. 
Recent studies do not support this claim. 
Massi et al. (1997) have shown that most 
molecular gas lies behind the Hll region 
and the cluster location is consistent with 
a face-on blister-type H ll region morphol- 
ogy. Optical spectroscopic study by Bohi- 
gas et al. (2004) indicates that the hot O 
stars of Pismis 24 emit 10^° UV photons 
s~^, sufficient to ionize the Hll region. 

Massey et al. (2001) consider the W- 
R/0 binary WR 93, lying 4' away from 
the cluster core, to be a likely member 
of Pismis 24 and caution that the clus- 
ter might be larger than the central con- 
centration of O stars. Spectroscopic study 
indicates a mass loss rate for WR 93 of 
M = 2.5 X 10~^ M0 yr~^ using a clumped 
wind model with terminal wind velocity 
= 2290 km s^^ (Prinja et al. 1990; 
Nugis et al. 1998; van der Hucht 2001). 
The nearby radio continuum and CO peaks 
suggest that WR 93 ionizes and heats the 
surrounding cloud (McBrcen ct al. 1983). 
The visual extinction towards WR 93 is 
high, with E{B - V) ^ 1.82 or Ay ~ 5.8 
mag (Nugis et al. 1998; van der Hucht 
2001). 

Healy et al. (2004b) and Hester & Desch 
(2005) present an HST/WFPC2 image of 
the G353.2+0.9 Hii region in NGC 6357, 
illustrating a sequence of formation of low- 
mass stars and describing their evolution 
from molecular cores to EGGs to proplyds 
in the vicinity of massive stars. Very re- 
cently, De Marco et al. (2006) analyze HST 



observations of eight HII regions, includ- 
ing Pismis 24, to look for protoplanetary 
disks. In Pismis 24, they find a jet-like 
feature, similar to a Herbig-Haro object in 
their source # 14. However, no central star 
is seen in the jet candidate. 

3. Observations and Data Reduc- 
tion 

3.1. Chandra Observation and Data 
Selection 

NGC 6357 was observed on July 9, 2004 
with the Imaging Array of the Advanced 
CCD Imaging Spectrometer (ACIS-I) on 
board Chandra. Detailed description of 
the instrument can be found in Weisskopf 
et al. (2002). Four front illuminated (FI) 
CCDs form the ACIS-I which covers a 
field-of-view (FOV) of ~ 17' x 17'. Two 
CCD chips on the ACIS Spectroscopic Ar- 
ray (ACIS-S) were also set to be func- 
tioning during the observation, although 
the mirror point spread function (PSF) 
degrades significantly far (> 15') off-axis. 
Here the ACIS-S data will not be quan- 
titatively discussed. The observation was 
made in the standard Timed Exposure, 
Very Faint mode, with 3.2-second integra- 
tion time and 5 pixel x 5 pixel event is- 
lands. The total exposure time is 38 ks 
and the satellite roll angle is 289 degrees. 
The aim point is centered on the 03 If star 
Pis 24-1, the heart of OB association Pis- 
mis 24. The observation ID is 4477. 

Appendix B of Townsley et al. (2003) 
and Getman et al. (2005a) provide de- 
tailed explanations of the Penn State ACIS 
data reduction procedure. Data reduction 
starts with filtering the Level 1 event list 
processed by the Chandra X-ray Center 
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pipeline to recover an improved Level 2 
event list. To improve absolute astrome- 
try, X-ray positions of AClS-1 sources were 
obtained by running the wavdetect wavelet- 
based source detection algorithm (Freeman 
et al. 2002) within the Chandra Interactive 
Analysis of Observations (CIAO) package 
on the original Level 2 event list, using only 
the central 8' x 8' of the field. The resulting 
X-ray sources were matched to the 2MASS 
point source catalog. We calculated the 
position offsets between 277 X-ray sources 
and their NIR counterparts, and applied 
an offset of -|-0".02 in right ascension (RA) 
and — 0".33 in declination (Dec) to the X- 
ray coordinates. The CCD charge transfer 
inefficiency (CTI) was corrected using the 
Penn State CTI corrector (version 1.16) 
developed by Townsley et al. (2000, 2002). 

3.2. Image Reconstruction and Source 
Finding 

A multifaceted source finding procedure 
is summarized here, which is designed to 
locate all potential sources, even in the 
presence of crowding or a broad PSF, with 
the recognition that features due to detec- 
tor noise may also enter the source list. 
First, we assemble a large number of can- 
didate sources using a variety of techniques 
and criteria (wavdetect, reconstruction, vi- 
sual selection, etc., as described below), 
which includes a number of possible false 
detections. These spurious sources are 
identified and removed after the photon ex- 
traction step (§3.3), where our customized 
tool evaluates the statistical significance of 
each detection above background, rejecting 
likely false detections when significances 
fail to meet a chosen threshold. These pho- 
tons are then returned to the background 



and we re-evaluate the source significance 
for the surviving candidates, finally defin- 
ing reliable detections and marginal detec- 
tions for inclusion in tables of source prop- 
erties. The details of this procedure are 
described below. 

Source detection was carried out with 
the improved Level 2 event list that is 
band-limited to 0.5-7.0 keV with cosmic 
ray afterglows removed. Effective expo- 
sure maps and images for the full field, the 
17' X 17' I array, the central 8' of the I ar- 
ray, and the inner 4' of the field were made 
with 4-pixel, 2.6-pixel, 1-pixcl, and half- 
pixel binning, respectively (an ACIS sky 
pixel is 0.5" on a side). Figure 2a shows 
the raw ACIS I-array image at reduced 
resolution. Adaptive-kernel smoothed fiux 
images were created with the CIAO tool 
csmooth (Ebeling et al. 2006) to help iden- 
tify additional potential faint sources and 
any possible diffuse emission. The result- 
ing full field two-color smoothed image and 
a zoom-in on the cluster core are shown 
in Figure 2; red and blue represent the 
soft band emission (0.5-2.0 keV) and hard 
band emission (2.0-7.0 keV), respectively. 

The wavdetect program was run with 
source significance threshold 1 x 10^^ 
on each of the four binned images de- 
scribed above, in three energy bands (0.5- 
2.0 keV, 2.0-7.0 keV, and 0.5-7.0 keV). 
This source significance threshold may give 
some false detections, but experience has 
shown us that this threshold is appropri- 
ate for crowded fields that might contain 
complex backgrounds due to diffuse emis- 
sion. Merging the detections from the 12 
runs resulted in a list of 614 potential point 
sources on ACIS-I. Careful visual inspec- 
tion revealed that wavdetect still missed 
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several apparent sources in the central re- 
gion of the cluster. To take advantage of 
the sub-arcsecond PSF at positions around 
the aimpoint, we applied a subpixel posi- 
tioning code (Mori et al. 2001) to improve 
spatial resolution in the inner part of the 
field and performed an image reconstruc- 
tion with the Lucy-Richardson maximum 
likelihood algorithm (Lucy 1974). 

Our image reconstruction technique for 
crowded stellar clusters is described in de- 
tail in Townslcy et al. (2006a) and a similar 
procedure is followed here. The maximum 
likelihood image reconstruction code ^ was 
run in the crowded central region using a 
~ 50" X 50" image made with 0.2-pixel bin- 
ning and centered on Pis 24-1, using the 
1.0 keV PSF at that on-axis location from 
the PSF library in the CI AO calibration 
database. We chose the 200th iteration 
of the algorithm as the most appropriate 
image for source searching, as it showed 
good PSF removal without over-resolving 
the data. We adopted the find procedure 
in the £)^OPi/Or package (Stetson 1987) 
to obtain the centroid, shape, and bright- 
ness parameters of the resolved peaks. Sev- 
eral combinations of flux (/), roundness 
(r), and sharpness (s) limits were tested 
to minimize the number of peaks from one 
photon, and the adopted parameters to 
reject spurious peaks were: r > 1.0 or 
r < —1.0; or / < 0.8; or / < 1.6 and s > 
0.98. Using SAOImage DS9 (Joye & Man- 
del 2003), world coordinate system (WCS) 
region files were made, centered on the un- 
rejected peaks, based on 90% PSF con- 
tours. The final reconstructed source list 



^Prom the IDL Astronomy User's Li- 
brary maintained by Wayne Landsman at 
http: / /idlastro. gsfc.nasa.gov/homepage. html. 



was made with the aid of the I- array im- 
age overlaid with these regions where any 
remaining reconstructed peaks that con- 
tain only single photon events within 90% 
PSF contours were rejected. This image re- 
construction procedure was performed sep- 
arately using another image centered on 
Pis 24-17 because Chandra's PSF varies at 
different locations and an appropriate lo- 
cal PSF must be used for the reconstruc- 
tion. Examples of valid sources from im- 
age reconstruction are shown in Figure 3. 
The list of potential sources from the maxi- 
mum likelihood reconstruction was merged 
with the wavdetect source catalog. By 
overlaying regions representing the poten- 
tial sources onto the original and smoothed 
ACIS-I images, an additional 40 sources 
were added based on spatial concentrations 
of > 3 photons and proximity of 2MASS 
near-IR sources (< 1"). The source find- 
ing procedure described here results in a 
total of 910 potential sources identified on 
the ACIS-I array. 

3.3. Photon Extraction and Limit- 
ing Sensitivity 

A preliminary event extraction for the 
910 potential X-ray sources was made with 
ACIS Extract (version 3.79; hereafter AE), 
a versatile script written in IDL that per- 
forms source extraction, fits X-ray spec- 
tra, creates light curves, and computes 
a wide variety of statistical properties of 
the sources (Broos et al. 2002)^. Using 
the ^£^calculated probability Pb that the 
extracted events are solely due to Pois- 
son fluctuations in the local background, 
source validity can be statistically eval- 

^http://www. astro. psu.edu/xray/docs/TARA/ae_users_guide. html 
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uated while taking account of the large 
distorted PSF at far off-axis locations 
and spatial variations in the background. 
The traditional source significance, de- 
fined as the photometric signal to noise 
ratio, is calculated for every source as 
well. We rejected sources with Pb > 1% 
likelihood of being a background fluctu- 
ation. The trimmed source list includes 
779 sources, with full-band (0.5-8.0 keV) 
net (background-subtracted) counts rang- 
ing from 1.7 to 1837 counts. We performed 
source extraction with this list using AE 
as described in detail by Townsley et al. 
(2006a). 

The 779 valid sources are purposely di- 
vided into two lists: the 665 sources with 
Pb < 0.1% make up our primary source 
list of highly reliable sources (Table 1), and 
the remaining 114 sources with Pb > 0.1% 
likelihood of being spurious background 
fluctuations^ are listed as tentative sources 
in Table 2. We believe that most of those 
tentative sources are likely real detections^. 

few Pb values arc greater than 0.01, our re- 
jection threshold in the preliminary run. This is 
due to the feedback to the background from the 
rejected sources when \vc re-extract the trimmed 
source list. As described in detail in §2.1 of Towns- 
ley et al. (2006a). our catalog cutting scheme is 
iterative: the initial set of 910 potential sources 
is extracted and then cut at the threshold of 1%. 
The rejected sources are allowed to contribute to 
the background and the trimmed list of 779 valid 
sources is then re-extracted, resulting in different 
local background estimates and different source 
validity metrics Pb- Thus, some Pb values in Ta- 
ble 2 are > 1%. 
^The fraction of tentative sources (Table 2) that 
have optical and near-infrared counterparts is 
50%, suggesting that they are not random spu- 
rious detections (§4.1). For the 64 weakest X-ray 
sources that have 3 counts or less, 10 of them have 
been cautioned with a "U" flag (unreliable sources 



We should note that the separation of 
reliable sources, tentative sources, and in- 
valid sources is based on a single simple 
statistic, the Poisson probability of ob- 
serving the extracted counts given the lo- 
cal background density. This permits sci- 
entists to evaluate the reliability of each 
source individually. A high cutoff (around 
4 counts for a typical Chandra expo- 
sure on-axis) reduces spurious noise-based 
sources {false positives) but misses some 
real sources {false negatives). A low cutoff 
reduces the false negatives while increasing 
the false positives. In the absence of a con- 
sensus statistical procedure for the Poisson 
weak-signal-in-background problem, there 
is no way to clean the list in an objec- 
tive fashion without losing real sources. 
Since different scientists make different de- 
cisions on the balance between false pos- 
itives and false negatives, and there are 
good scientific reasons to believe many as- 
tronomically interesting faint sources lie 
in Chandra images of rich stellar clusters, 
we have opted to report marginal sources 
in Table 2. The reader is clearly warned 
that the reality of some of these weak- 
est sources is suspect, and each source's 
Pb can be examined individually. These 
choices concerning how we choose to report 
faint sources do not affect our science re- 
sults. The weakest X-ray sources are only 
used for X-ray luminosity function estima- 
tion (§5.1), contributing to the faintest two 
bins of the XLF, which do not impact on 
our conclusions regarding the NGC 6357 

that have large probability of being background, 
§4.5) in the table, while 27 of them have optical 
and near-infrared identifications. Most of these 
faint sources are clustered within the central 2' of 
the field and are likely low-mass cluster members. 
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stellar population. We further note that 

every detection scheme must - explicitly or 
implicitly - estimate the local background, 
and every such estimate must implicitly 
identify the data that are guessed not to 
be background. For crowded young stellar 
clusters, this can never be accomplished 
in a unique, objective manner as it de- 
pends on the unknown distribution of true 
faint sources (including thousands of stars 
which lie below any detection threshold) 
and possibly diffuse emission from O star 
winds (Townsley et al. 2003). 

Table 1 and Table 2 have a format that 
closely follows the tables in Getman et al. 
(2006) and are identical to Tables 1 and 
2 in Townsley et al. (2006a). Column 1 
gives the ACIS running sequence number, 
and column 2 provides the lAU designa- 
tion. Columns 3-6 show RA and Dec in 
degrees, positional error in arcseconds, and 
off-axis angle in arcminutes. Column 7- 
11 give net counts in the full band, the 
associated Ici-equivalent uncertainty based 
on Poisson statistics, background counts in 
the extraction region, net counts in the 
hard band, and the fraction of the PSF 
used for source extraction (reduced val- 
ues for this quantity indicate source crowd- 
ing). Columns 12 and 13 present two mea- 
sures of the source validity: the photomet- 
ric significance, and the probability that 
the source is a spurious background fluc- 
tuation Pb, used for establishing source 
validity as described above. Column 14 
lists possible anomalies in the observation 
and Column 15 checks the source variabil- 
ity. Table footnotes give details for Col- 
umn 13-15. Column 16 gives the "effec- 
tive" exposure time, defined as the amount 
of exposure time needed for the source to 



accumulate the current number of counts 
if placed at the Chandra aimpoint. Col- 
umn 17 gives the source's median energy 
in the full spectral band after background 
subtraction. 

The on-axis limiting sensitivity of this 

38 ks observation can be estimated us- 
ing the Portable Interactive Multi-Mission 
Simulator (PIMMS ^). Assuming an on- 
axis detection of 3 counts (0.5-8.0 keV) 
in 38170 seconds and an absorbing column 
Nh — 1.0 X 10^^ cm~^ (estimated from 
an average Ay ~ 6.0 mag to Pismis 24 
in Bohigas et al. 2004, adopting the con- 
version between Nh and Ay from Vuong 
et al. 2003), the observed X-ray flux for 
a thermal plasma with temperature of 1 
keV (appropriate for a pre-MS star's X- 
ray spectrum) is 5.6 x 10^^^ ergs s^^ cm~^. 
Correcting for absorption and a distance of 
2.56 kpc to Pismis 24, this gives a limiting 
luminosity of logL^ = 30.2 erg s~^, where 
Lf c is the intrinsic (absorption-corrected) 
luminosity in the total Chandra band 0.5- 
8.0 keV. 

3.4. Source Vciriability and Spectral 
Fitting 

A band-limited, adaptively smoothed 
light curve for each source was generated 
by yli? together with a median energy time- 
series. The variability of the source was 
evaluated by the significance of a one-sided 
Kolmogorov-Smirnov statistic Pks, com- 
paring the source events arrival times to 
that of a uniform light curve model. If 
Pks > 0.05, the null hypothesis of a con- 
stant source was accepted. If Pks < 
0.005, the source was classified as "vari- 

^http://heas arc. gsfc.nasa.gov/Tools/w3pimms. html 
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able." Otherwise, the source was "possibly 
variable." 

For the 40 sources with >80 net counts, 
the extracted spectra were fit to models of 
optically thin thermal plasmas using the 
XSPEC pachsige (Arnaud 1996), based on 
source spectra, background spectra, an- 
cillary response functions (ARFs) and re- 
distribution matrix functions (RMFs) con- 
structed by the AE script. The data were 
grouped in bins with > 10 counts and the 
statistic was adopted to evaluate the 
goodness-of-fit. A single apec (Smith et al. 
2001) thermal plasma model with a sin- 
gle absorption component was considered. 
Interactive fitting with alternative models 
(e.g., vapec thermal plasma with variable 
abundances; two-temperature apec ther- 
mal plasma model) was performed for a 
few strong sources that did not achieve a 
satisfying fit with the automated process. 
Table 3 gives the best fit results using the 
thermal plasma model; its format closely 
follows Table 3 in Townsley et al. (2006a). 
Uncertainties representing 90% confidence 
intervals are given; incomplete confidence 
levels imply the parameters are not well 
constrained, or XSPEC may have encoun- 
tered some abnormality in the error cal- 
culations. For the 144 sources that have 
net counts between 20 and 80, we fit un- 
grouped spectra with a single temperature 
apec thermal plasma model using the like- 
lihood ratio test (C statistic); the fitting 
results are presented in Table 4. 

None of the NGC 6357 sources are 
bright enough to warrant correction for 
photon pile-up. The brightest source in the 
field is the W-R/0 binary WR 93 which 
has 1837 ACIS counts in the 38 ks expo- 
sure, or a count rate of 0.15 counts per 



CCD frame. As it hes 5'. 4 off-axis where 

the PSF is somewhat degraded, the photon 
pile-up effect will not be significant enough 
to distort the spectral fitting, according to 
previous simulation work (Townsley et al. 
2002). 

3.5. SIRIUS NIR Observation 

To identify infrared counterparts with 
images that are deeper and of higher spa- 
tial resolution than 2MASS, J, H, and K 
band images were obtained on 30 Septem- 
ber 2005 using the wide-field NIR SIR- 
IUS camera on the 1.4 m IRSF telescope 
at the South Africa Astronomical Obser- 
vatory. Details of the instrument can be 
found in Nagayama et al. (2003). The pixel 
scale of the detector is 0.45 arcsec pixel"''" 
with a 7'. 7 x 7'. 7 field-of-view. Three sec- 
ond exposures were made at each of 20 
dither positions for a total exposure of 60 
seconds. A typical limiting magnitude for 
such SIRIUS observations is ^^=15. 7. Un- 
fortunately, the observing conditions were 
non-photometric and the seeing was 1.1" 
FWHM in the K band. We nonetheless 
found these images to be very valuable for 
identifying infrared counterparts (as de- 
scribed below) because of the higher spa- 
tial resolution compared to 2MASS. 

4. Identification of the Chandra 
Sources 

4.1. X-ray Sources with Stellar Coun- 
terparts 

Identifying stellar counterparts for X- 
ray sources in massive star forming regions 
beyond 1 kpc is challenging. In regions 
such as the Chamaeleon I cloud (0.2 kpc) 
or Cepheus B region (0.75 kpc), the cluster 
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members are sufficiently bright and well- 
separated that the 2MASS catalog pro- 
vides adequate identifiers for nearly all 
Chandra stellar sources (Feigelson & Law- 
son 2004; Getman et al. 2006). However, 
for NGC 6357 at 2.5 kpc, deeper and higher 
resolution observations such as the SIRIUS 
images are needed for source identification. 
In massive star forming regions such as 
NGC 6357, optical and NIR imaging have 
the additional challenge of detecting faint 
stars amidst highly structured background 
emission from the H II region or deeply em- 
bedded in the molecular cloud. Thus we 
cannot arrive at a quantitative statement 
regarding the completeness of the search 
for stellar counterparts of Chandra sources. 
But we can say that, with very few excep- 
tions, e.g., infrared luminous active galac- 
tic nuclei (AGNs), the counterparts which 
we do identify are indeed cluster members. 

For sources with off-axis angle 9 < 3', 
the Chandra sourcelist was matched to the 
2MASS Catalog of Point Sources with a 
1" matching radius. The search radius 
was enlarged to 2" for ^ > 3' in con- 
sideration of the degraded off-axis Chan- 
dra PSF (Cetman et al. 2005a). We find 
that 445 of our X-ray sources (58%) have 
2MASS counterparts. The matching re- 
sults are presented in Table 5. Four of the 
ACIS sources (#88, 517, 688, 774) have 
two counterparts within 2" and we assign 
the star with the smaller separation as the 
real counterpart. Three pairs of sources - 
252 and 257, 401 and 407, 404 and 406 - 
share a single 2MASS counterpart within 
1". Given the high density of both X-ray 
sources and 2MASS objects around cluster 
center, we follow the simple Monte Carlo 
technique from Alexander et al. (2003) to 



evaluate the number of chance superposi- 
tions between 2MASS sources and X-ray 
detections. A random offset between 5" 
and 10" is applied to the catalog to remove 
the true associations, and then the same 
cross-identification radii are used to match 
X-ray positions to 2MASS positions. This 
procedure is repeated many times and the 
resulting false detection rate is 2 — 3%. 

For the 564 X-ray sources covered by 
SIRIUS imaging, we visually identified 476 
infrared counterparts with separations < 
1" from X-ray sources. Of these, 308 
are known 2MASS counterparts that have 
been matched to ACIS sources. Therefore 
the combined matches from 2MASS and 
SIRIUS yield 613 NIR counterparts for the 
779 X-ray point sources. Of the 166 X-ray 
sources without NIR counterparts, 89 lie in 
the SIRIUS FOV. 

In spite of its small FOV, an archival 
HST WFPC2 (F814W) image of the Hll 
region G353.2+0.9 in NGC 6357 (Figure 4; 
Healy et al. 2004b) was visually examined 
for optical counterparts. By registering to 
a 2MASS image, an astrometric correction 
of ARA = 0.9" and ADec = -0.4" was ap- 
plied to the HST image. There are 154 X- 
ray sources covered by the WFPC2 imag- 
ing, 99 of which have HST counterparts. 
The WFPC2 image demonstrates details of 
this interface between the hot tenuous H ll 
region and the cold structured molecular 
cloud. Some of the optically bright sources 
are not X-ray emitters and are probably 
foreground stars; a few stars appear very 
faint and are probably deeply embedded in 
the molecular cloud. All but three of the 
HST counterparts (corresponding to ACIS 
#359, 375, 390) have matched infrared 
counterparts. This increases the number 
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of X-ray sources with ONIR counterparts 
to 616. Thanks to HST's high resolution 
imaging, those ACTS sources are resolved 
from the crowded field where infrared im- 
ages fail to detect them. Even though opti- 
cal imaging suffers heavily from extinction, 
we speculate that ~ 10 — 20 more HST 
counterparts for the less-embedded ACIS 
sources with separations < 1" from other 
sources would have been recovered due to 
the high spatial resolution, if the HST spa- 
tial coverage of Pismis 24 were complete. 

4.2. Extragalactic Contamination 

We estimate the contamination from 
the extragalactic population following the 
method described in Getman et al. (2006). 
Monte-Carlo simulations of the expected 
extragalactic source contamination are 
constructed by placing fake AGNs ran- 
domly across the ACIS-I field superposed 
on the AGIS background. The X-ray flux 
distribution follows the logA^-logS distri- 
bution from Moretti et al. (2003) and the 
corresponding power law photon index for 
individual sources follows the flux depen- 
dence found by Brandt et al. (2001). Syn- 
thetic source spectra were generated with 
the XSPEC fakeit function and convolved 
with an absorption column Nh- We ran 
simulations with three typical Nh- the to- 
tal Galactic absorption towards NGG 6357 
{Nh = 1.2 X 10^2 cm-2. Dickey & Lock- 
man 1990); Galactic absorption plus an 
additional Nh = 1-0 x 10^^ cm"^ typical of 
Pismis 24 OB stars; and Galactic absorp- 
tion plus an additional Nh — 3.0 x 10^^ 
cm~2 typical of reddened 2MASS stars. 
Then we identified the closest real X-ray 
detection to the fake source. Source sig- 
nificance for each fake AGN was then cal- 



culated using the local background of the 
nearby real source in the NGC 6357 back- 
ground map. We set a detection threshold 
of Signif > 2.3, where Signif is the quan- 
tity given in column 12 of Table 1. This 
limiting source significance is the median 
Signif of those hard X-ray sources with- 
out infrared counterparts and lying far off- 
axis, which are the most likely candidates 
for AGNs in the field. 

The number of detectable AGNs from 
typical simulations with the three absorp- 
tion column densities ranges from 14 to 
29. The level of AGN contamination for 
the NGC 6357 region (Galactic coordinates 
I — 353.2, b — +0.9) is similar to the con- 
tamination calculated for the Cepheus B 
(/ = 110.2, b = +2.6) observation which 
had a similar Chandra exposure (Getman 
et al. 2006). The true extragalactic con- 
tamination may be lower due to local- 
ized regions of very high absorption by 
the molecular cloud behind the OB asso- 
ciation. We conclude that no more than 
4%, and more hkely ~ 2%, of the 779 
X-ray sources are extragalactic contami- 
nants. Among possible AGN candidates, 
the brightest in X-rays is source ^1. Based 
on its fiat X-ray spectrum, large separa- 
tion from the central cluster, non-variable 
hghtcurve, and lack of optical or near-IR 
counterpart, AGIS #1 is considered to be a 
likely AGN. Using a power law model with 
a single absorption component, its spec- 
trum can be adequately fit using a power 
law photon index F = 1.4 and an absorp- 
tion column Nh — 8.0 x 10^^ cm"^. 

4.3. Field Star Contamination 

To estimate possible contamination 
from foreground field stars, we follow 
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the techniques described in the Chandra 
Orion Ultradeep Project (COUP) mem- 
bership study and the Chandra study of 
the Cep B star forming region (Getman 
et ah 2005b, 2006). We use simulations 
based on the stellar population synthesis 
model of Galactic disk stars by the Be- 
sancon group (Robin ct al. 2003). The 
models are calculated using their on-line 
service^. Synthetic catalogues of stars gen- 
erated from these models predict ~ 5600 
field stars with V < 22 mag in a solid an- 
gle of ~ 0.08 deg2 (AGlS-1 FOV) to the 
distance of NGG 6357. About 10% are F- 
type main sequence (MS) stars, ~ 20% are 
G-type MS stars, ~ 45% are K-type MS 
stars, ~ 20% are M-type dwarfs; the rest 
are giants with a few subdwarfs. 

We then adopt the X-ray luminos- 
ity functions (XLFs) for MS K-M stars 
(Schmitt et al. 1995), F-G stars (Schmitt 
1997), and late type giants (Huensch et 
al. 1996) derived from a complete volume- 
limited sample in the solar neighborhood. 
The conversion of intrinsic X-ray flux in 
the ROSAT PSPG 0.1 - 2.4 keV band 
to Chandra AGIS-I 0.5 - 8.0 keV count 
rate is from PIMMS, using a thermal 
plasma model assuming kT ^ {And'^Fx x 
5.5 X 10"26)0-2/11.34 keV from Giidel et 
al. (1998), where d is the stellar distance 
and Fx is the X-ray flux in ROSAT PSPG 
0.1—2.4 keV band. Only Galactic absorp- 
tion is assumed to be present for stars in 
the foreground of the NGG 6357 cloud and 
the amount of absorption to the fake star 
is calculated based on its distance when 
the simulated catalog is generated. 

Following the same procedures used in 

^http://bison.obs-besancon.fr/niodele/ 



the simulations for extragalactic contami- 
nation, a simulated field of ~ 5600 X-ray 
field stars was added to the ACIS back- 
ground map, and the star locations were 
examined for source significance. For de- 
tection thresholds ranging from Signif ~ 
1.5 to Signif ~ 3.0, 10000 simulation 
runs give a typical number of 6 to 16 
detectable foreground field stars. These 
sources should be readily distinguished 
from extragalactic AGN and most young 
cluster members by the brightness of their 
stellar counterparts and (when sufficient 
X-ray counts are present to measure me- 
dian energies) low absorption in their spec- 
tra. Eleven individual stars are identified 
as likely foregound stars in the field (§4.5). 
Those contaminants will be removed from 
further consideration of the young cluster 
population. The hkelihood of background 
star contamination is small because of the 
heavy absorption of the molecular cloud 
and the large distance. 

The contamination from Galactic disk 
stars is thus 1 — 2% of the observed source 
population, about half that of the extra- 
galactic contamination. A complete re- 
moval of all non-members is impossible 
without further information, thus we ig- 
nore the remaining foreground star and 
AGN contamination in our further anal- 
ysis given that the fraction is statistically 
very small. 

4.4. Likely New Stellar Members 

For the 163 ACIS sources without ONIR 
counterparts, at most 20-30 of them are 
possible extragalactic contaminants (§4.2). 
The rest of the ~140 sources are proba- 
bly new members of Pismis 24, in addition 
to the 605 likely members (excluding the 
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11 candidate foreground stars) with ONIR 
counterparts identified from the thousands 
of ONIR stars in the field (§4.1). 

Sixty-two ACIS sources without ONIR 
counterparts are clustered within 2' of the 
03 star Pis 24-1. The rest are widely dis- 
tributed throughout the field, although a 
few are located to the northwest of the ion- 
ization front where there is much less opti- 
cal nebulosity. Five of them are clustered 
around the WR star. Most of these likely 
members are faint X-ray sources; only two 
of them have more than 10 net counts. 
More than 100 of them have median pho- 
ton energy {E) > 2.0 keV and ~ 50 have 
{E) > 3.0 keV, implying heavy obscura- 
tion (see Figure 8 from Feigelson et al. 
2005). We speculate that, due to the face- 
on geometric configuration, these repre- 
sent the more distant cluster members that 
are located behind those optically visible 
members, deeply immersed in the molecu- 
lar cloud. Some may be very young, em- 
bedded in molecular cores or protostellar 
envelopes that have not yet been identi- 
fied. But many are likely to be older Class 
II and III pre-MS stars that have counter- 
parts in Spitzer Space Telescope surveys. 
The properties and spatial distribution of 
these likely new members will be further 
discussed in §5.2 and §6.5, after we classify 
the X-ray sources. 

4.5. Classification of tiie X-ray Sam- 
ple 

From the results of the previous sec- 
tions, we can establish a rough classifica- 
tion of the X-ray sources. These class des- 
ignations appear in Table 5. The numbers 
in parentheses after the class names are the 
number of sources in each category. 



R (600)- Rehable X-ray sources with 
ONIR stellar identifications. They 
represent the more massive and less 
obscured members of the stellar pop- 
ulation. 

C (590)— X-ray sources with off-axis an- 
gle < 5.0'. They are likely mem- 
bers of the central cluster. 

H (109)— Hard X-ray sources with {E) > 
3.0 keV. Most of these are heav- 
ily obscured members of the embed- 
ded stellar population. About 20-30 
of these are probably extragalactic 
AGNs. 

U (18)- Unreliable X-ray detections. We 
consider the 18 X-ray sources in Ta- 
ble 2 with large probability of being 
background (logPg > —1.8) to be 
possible spurious detections. 

F (11)— Foreground field stars. Seven 
soft X-ray sources {{E) < 1.3 keV, 
corresponding to log Nh < 22.0 
cm~^) with more than 20 net counts 
in 0.5-8.0 keV (to have reliable me- 
dian photon energy measurements) 
are identified as foreground sources, 
excluding any known O stars. One 
X-ray emitting A8 subgiant, HD 
157528, was recovered using this cri- 
terion. The limitation for this cri- 
terion is that unidentified massive 
stars with soft X-ray emission may 
be included. Two of these candidate 
foreground stars, ACIS #9, #539, 
are included as candiate O stars (see 
§6.2). In addition, we select 4 soft 
X-ray sources {{E) < 1.3 keV) as 
likely field stars since their JHK col- 
ors are consistent with the locus of 



14 



Besangon simulated field stars. All 

these sources have visual photometry 
in the UCAC2 catalogue (Zacharias 
et al. 2004), which provides addi- 
tional support for them being classi- 
fied as hkely foreground stars. 

The classes R, C, and H are not ex- 
clusive and many sources fall in two or 
three categories as illustrated in Figure 5. 
Thirty-eight sources are not classified be- 
cause they are either too faint or poorly 
characterized. 

Moffat & Vogt (1973) listed 15 bright 
stars in Pismis 24 as their photometry tar- 
gets and concluded 12 of them are mem- 
bers. Spectroscopic observation confirmed 
that seven of these, and four additional 
stars, are massive OB stars (Massey et al. 
2001). Our X-ray sample with infrared 
counterparts increases the known cluster 
population 50-fold, from ~ 15 to ~ 750, 
providing the first opportunity to investi- 
gate the full stellar population of the re- 
gion. 

The H—Kws. J—H color-color diagram 
of 142 Chandra sources with high-quality 
2MASS photometry is shown in Figure 6. 
The typical uncertainty in color is ~ 0.06 
mag. The locus for main-sequence stars 
of different spectral types from early O to 
late M (Bessell & Brett 1988) is plotted 
as the thick fine. Stars that are reddened 
by normal interstellar extinction occupy 

■^The 2MASS sources plotted here and in the J - 
vs. J diagram have signal-to-noise ratios greater 
than 6 (photometric quahty higher than "C" ) and 
are free of confusion flags and warning of contam- 
ination flags. Sec Explanatory Supplement to the 
2MASS AU Sky Data Release (Cutri et al. 2003; 
http: / /www.ipac.caltech.edu/2mass/releases/allsky/d 
for details. 



the color space between the two redden- 
ing vectors of an O dwarf and an M giant 
(dashed lines). The foreground AS sub- 
giant (# 773) clearly reveals itself, lying on 
the locus with almost no reddening. The 
stars to the right of the reddening band 
are 7^-band excess sources. It is clearly 
seen that the X-ray observation selects far 
more Class 111 infrared sources (stars that 
have simple blackbody spectral-energy dis- 
tributions, with little or no inner accretion 
disks) than Class II sources (stars with K- 
band excess circumstellar disks; Lada 1987; 
Wolk & Walter 1996). 

The 2MASS J—H vs. J color-magnitude 
diagram (Figure 7) gives both the rough 
mass distribution and the absorption dis- 
tribution of the ACIS sources with 2MASS 
counterparts. Zero age main sequence 
(ZAMS) track (Cox 2000) and two pre- 
MS isochrones (Baraffe et al. 1998; Siess 
et al. 2000) are shown with spectral types, 
stellar masses, and reddening vectors for 
estimation of the dcreddened masses and 
the associated absorption. We assume 1 
Myr for the age of Pis 24 cluster and a dis- 
tance of 2.56 kpc when estimating stellar 
properties from the color magnitude dia- 
gram. We remind the reader that the esti- 
mated mass and the amount of absorption 
for a given star also relies on the assump- 
tion of distance and cluster age besides 
the uncertainties in the star's color and 
magnitude. For example, if we adopt a 
younger cluster age (0.3 Myr), the pre-MS 
isochrone will shift upwards and yield a 
lower dereddened mass and less amount of 
absorption as shown in Figure 7. Note that 
the inferred properties for stars with spec- 
tral type later than B3 are less dependent 
oc/)on the age and distance assumptions. In 
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addition we show the H — K vs. K color- 
magnitude diagram for the same sample 
of stars accounting for some small i^'-band 
excess. Known Pis 24 optical members de- 
tected in X-ray are located along the red- 
dening vectors from the locus of O stars 
with a typical Ay ~ 5 mag. Not surpris- 
ingly, the foreground star 7^ 773 resembles 
a very massive and luminous object at the 
distance of NGC 6357. Besides the fore- 
ground star, WR 93 is the most luminous 
object. 

Most of the other ACIS sources occupy 
infrared color domains that arc consistent 
with being high-mass stars and pre-MS 
intermediate-mass stars (M > 2Mq) with 
5 < Ay < 15 at 2.5 kpc. The highly- 
reddened bright 2MASS sources are pos- 
sibly obscured massive stars. The new 
high-mass stars double the known high- 
mass population of the cluster (§6.2). The 
intermediate-mass stars, mostly Class III 
objects, are newly identified members of 
the cluster (§6.3). The ratio between Class 
111 and Class II objects for the high-mass 
and intermediate-mass cluster members is 
~ 15 : 1 based on the ~140 stars with 
mass M > 2Mq. Although the uncertain- 
ties of the estimated stellar masses are not 
negligible (as much as a factor of 2; see 
the discussion of Orion stars in Preibisch 
et al. 2005; Luhman 1999; Hillenbrand & 
White 2004), this is a quantitative mea- 
surement for the massive end of the clus- 
ter members and implies that the fraction 
of stars possessing inner disks among the 
intermediate-mass and high-mass stars is 
low even in such a young stellar cluster. 
Those intermediate-mass stars with signif- 
icant IR excess are likely the accreting pre- 
MS stars known as Herbig Ae/Be stars 



(HAeBes; Herbig 1960; The et al. 1994). 

X-ray properties of the intermediate mass 
stars with infrared excess will be examined 
in § 6. A handful of objects are heav- 
ily obscured with Ay > 20. Due to their 
brightness in K-band (X^ ~ 9 — 13), these 
objects, some of which exhibit K-band ex- 
cess, are easily accessible for ground based 
follow-up infrared spectroscopic or photo- 
metric study. The remaining ~ 170 ACIS 
members without 2MASS counterparts but 
seen with SIRIUS and/or HST are likely 
lower mass stars. The close separations 
to the high/intermediate mass stars sug- 
gest they perhaps are companions, but this 
could well be a spatial projection effect in 
such a crowded cluster. Similar to those 
unobscured ACIS sources without ONIR 
counterparts, they are either too faint or 
too close to bright stars to be cataloged in 
2MASS. 

4.6. EGGs and Protostars 

Two-thirds of the IR counterparts - 
those with low photometric quality - do 
not appear in Figure 6. Although candi- 
date protostars may be present, the large 
uncertainty in their IR color and the upper 
limits in i^-band do not allow us to con- 
fidently classify class II/I objects. Precise 
infrared photometry will be the key to es- 
tablishing whether they are true embedded 
protostellar objects. 

When compared to an archival HST 
WFPC2 image (Figure 4; Healy et al. 
2004b), the X-ray source #461 (CXOU 
J172445. 74-341106.9) is spatially coinci- 
dent with the tip of an evaporating clump 
extending from the wall of the molecular 
cloud, which is known as an evaporating 
gaseous globule (EGG; Hester et al. 1996; 
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Hester & Desch 2005). It has only 4.6 

net counts but the X-ray photons are spa- 
tially concentrated; it is a reliable Table 1 
source. The X-ray emission appears to 
be moderately absorbed with {E) =2.5 
keV. The colors of its infrared counterpart 
show large NIR excess {J — H = 1.06, 
H - K = 1.42). The 2MASS counterpart 
was previously detected in the Persi et al. 
(1986) infrared observation (=IRS 4) and 
its JHKLMNQ spectral energy distribu- 
tion (SED) gives a 1-20 fim luminosity of 
630 Lq (Persi et al. 1986). This luminosity 
and its infrared colors are consistent with 
that of a protostar with spectral type BO- 
B2 (Persi et al. 1986; Bohigas et al. 2004). 
In addition, it also coincides with the com- 
pact radio peak of a 6 cm VLA observation 
(=VLA A; Felh et al. 1990). 

CXOU J172445. 74-341106.9 may be the 
first detection of X-ray emission from an 
EGG, which represents one important evo- 
lutionary stage of star formation. Linsky 
et al. (2003) identified deeply embedded 
hard X-ray sources in the "pillars of cre- 
ation" of the Eagle Nebula (M 16; Hester et 
al. 1996) but were unable to find any such 
sources associated with the EGGs. The 
non-detection of X-ray emission in M 16 
EGGs may be related to their high thresh- 
old of detection (>6 counts), heavy extinc- 
tion, or the spectral type and evolutionary 
phases of the embedded stars. Although 
the uncertainty in the X-ray luminosity 
is potentially very large for such a faint 
source, the estimated ~ 30.1 ergs s~^ 
from PIMMS is somewhat consistent with 
other X-ray detected high-mass protostars 
(NGC 2024 IRS 2b, 08-B2, Skinner et al. 
[2003]; S106 IRS4, 07-BO, Giardino et al. 
[2004]). It is likely that the actual X-ray 



source is the embedded protostar instead 

of the gaseous globule. 

No X-ray detection is associated with 
the HST jet candidate, Pismis 24 ^ 14 in 
De Marco et al. (2006). We examined the 
X-ray image at the HST position and no 
photons can be found within a 90% PSF 
contour. 

5. Properties of the Stellar Cluster 
and Its Environment 

5.1. X-ray Luminosity Function 

We construct the observed hard band 
(2-8 keV) XLF for all X-ray sources, ex- 
cluding sources classified as U and F. 
Based on the fact that the extragalactic 
contamination is small, we only exclude 
ACIS # 1 that is very likely an AGN. We 
do not attempt to study the observed total 
band (0.5-8.0 keV) XLF or absorption cor- 
rected total band luminosities {Lt^c: 0.5-8.0 
keV), since we do not have reliable mea- 
surements of Nh for faint sources. In the 
case of heavily absorbed X-ray sources, 
the unknown soft component can intro- 
duce a large uncertainty in the observed 
total band X-ray luminosity. 

For sources with more than 20 net 
counts, the hard band (2.0-8.0 keV) fiux 
Fh and absorption-corrected hard band 
fiux Fh^c (and therefore the luminosities 
Lfi and Lh^c) are available from XSPEC 
spectral fitting. For fainter sources, we 
define the following procedure to estimate 
Fh and also F/j.^ consistently. First, if the 
faint X-ray source has a 2MASS coun- 
terpart with reliable photometry, we use 
the Ay derived from the color magni- 
tude diagram to obtain the Nh, assum- 
ing a 1 Myr age for the cluster mem- 
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bers Nh/Av = 1.6 x 10^^ cm-^ mag~^ 
(Vuong et al. 2003). We note that adopt- 
ing a younger age or a smaller distance 
will shift the pre-MS isochrone upwards 
in the color magnitude diagram, causing 
generally lower values for the derived Ay- 
However for these faint sources that are 
likely low-mass cluster members, the de- 
crease in Ay is very small, ~ 1 mag (Fig- 
ure 7) . Otherwise, as shown in Feigelson et 
al. (2005), empirical conversion from me- 
dian photon energy {E) provides us with 
an effective approximate logNu value for 
sources with only a few counts: XogNn ~ 
9.96 + 13.62 X (E) - 3.86 x {E)^ cm'^ for 
relatively soft sources (1.0 < {E) < 1.7 
keV); log Nh = 21.22 + 0.44 x (E) cm'^ 
for harder ones {{E) > 1.7 keV). We esti- 
mate that the uncertainties in logiVff are 
~ ±0.3 and can be as large as ±0.8 for soft 
sources that have (E) < 1.3 keV. We as- 
sign plasma temperature kT based on {E) 
to reflect the observed trend of increasing 
temperature with increasing median en- 
ergy: typically kT — 1 keV for {E) < 1.3 
keV, kT ^ 2 keV for 1.3 < (E) < 2.0 
keV, kT = 3 keV for 2.0 < (E) < 4.0 keV, 
and kT = 5 keV for {E) > 4.0 keV (see 
Figure 8 in Feigelson et al. 2005). With 
the Nh and kT parameters set and ACIS-I 
count rates known for each source, we run 
PIMMS simulations to derive the observed 
(absorbed) and absorption corrected flux 
in the hard band. 

The XLFs for the lightly obscured 
{{E) < 3.0 keV) and heavily obscured 
{{E) > 3.0 keV) stellar populations in 
NGC 6357 are plotted in Figure 9. For 
comparison, XLFs derived for the lightly 
obscured COUP Orion Nebula population 
(Feigelson et al. 2005) and for the Cep 



B unobscured and obscured populations 

(Getman et al. 2006) are also shown. The 
shapes of the high luminosity portions of 
the lightly obscured and obscured COUP 
XLFs (see Figure 3 in Feigelson et al. 2005) 
are very similar, but the weaker sources in 
the obscured COUP sample are missing. 
Therefore we adopt the lightly obscured 
COUP XLF as the template in this popu- 
lation cahbration exercise. At the high Lx 
end, the NGC 6357 XLF is clearly consis- 
tent with the power law relation as seen 
in COUP and Cep B, allowing us to scale 
to the COUP XLF and make the first re- 
hable estimate of Pismis 24 's total cluster 
population. 

At logL;, ~ 30.3 ergs s-\ the NGC 6357 
XLF begins to flatten and decline dramati- 
cally. We examine the completeness in the 
source detection to determine whether the 
decline is due to a real deficiency of cluster 
members in this luminosity range or is sim- 
ply an expected effect when reaching the 
faint source detection limit. If the observed 
X-ray flux from a cluster member with 
such a luminosity is well above our sensitiv- 
ity limit, the deficiency of stars in this lu- 
minosity range would imply a different X- 
ray luminosity distribution for the popula- 
tion. Given that sources with 10 net counts 
in the full-band at any off-axis angle are 
confidently detected, we adopt two differ- 
ent methods to estimate what is the corre- 
sponding flux limit to 10 net counts. First, 
assuming a 2 keV thermal plasma, PIMMS 
simulation gives logF^ — —14.6 ergs s~^ 
cm~^ for a count rate of 10/38 counts ks~^. 
Independently, using all sources with more 
than 20 net counts that have flux estimates 
from spectral fitting, we derive an empiri- 
cal conversion factor between fiux and net 
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counts for this observation, and the 10 net 
counts also converts to log Fh ~ —14.6 
ergs s~^ cm~^. Therefore the two ap- 
proaches give the same Fh, which converts 
to logL^ ~ 30.2 ergs s~^. The dechne in 
the NGC 6357 XLF below logL^ ~ 30.3 
ergs s~^ is thus consistent with incomplete 
detection of the less luminous cluster mem- 
bers. We do not have sufficient sensitiv- 
ity to search for a difference in XLF shape 
compared to COUP below log ~ 30 ergs 
s~^, as seen in Cep B by Getman et al. 
(2006). 

The scaling factor to match the COUP 
XLF to the NGC 6357 XLF is ~ 5 - 6 
at higher values of log Lh for the lightly 
obscured population. Using the XLF con- 
structed with Lh,c gives the same scaling 
factor. The XLF from Lh for the obscured 
stellar population in NGC 6357 is compa- 
rable to the COUP XLF, while the ab- 
sorption corrected XLF seems to be con- 
sistently ~ 5 times higher than the COUP 
XLF for the high luminosity bins. Alto- 
gether, this implies that the X-ray emitting 
stellar population in NGC 6357 is roughly 
5 — 10 times^ richer than the lightly ob- 
scured Orion Nebula Cluster (839 unob- 
scured cool stars; Feigelson et al. 2005). 
Therefore the X-ray emitting population 
in Pismis 24 is ~ 5000. Given that the 
ONC has ~ 2000 known members from the 
deep optical study by Hillenbrand (1997) 
and COUP detects more than half of the 
stellar population in X-rays, the total stel- 
lar population in Pismis 24 is estimated to 

'If the true distance to the cluster is a smaller value, 
e.g. 1.7 kpc, all the derived luminosities, including 
logLfi and logL;i,c, will be systematically lower by 
0.35 dex and the estimated population will be ~3 
times the Orion Nebular Cluster. 



be two times larger than the population 

detected in X-rays, ~ 10, 000 members. It 
is not unexpected, assuming the standard 
IMF, to find massive stars with > lOOM© 
in a cluster with such a large population. A 
note of caution here is that, even in such a 
young star-forming region, the stellar pop- 
ulation in consideration is not perfectly co- 
eval. It perhaps includes a younger pop- 
ulation triggered to form by the ionizing 
front from the massive stars, but the small 
population of the second generation stars 
has little impact on the IMF of this rich 
cluster. 

5.2. Spatial Distribution of the Stel- 
lax Cluster 

Aside from several dozen known and 
new OB stars (see §6), the locations of over 
700 new lower-mass cluster members arc 
now determined with sub-arcsec precision. 
Figure 10 shows the spatial distribution 
of the cluster members. The large scale 
optical nebulosity contour from the DSS 
R-band image, which also traces the mid- 
IR ring-like morphology, is shown together 
with the outlines of the ionization front in 
yellow and the CO "South-Eastern Com- 
plex" in green (S.E.C. in the CO emission 
map in Massi et al. 1997). The geometry 
of the cluster appears spherically symmet- 
ric: the source density in the cluster core 
near the ACIS-I aimpoint is the highest, 
and declines radially with no apparent dis- 
continuity. There seem to be fewer X-ray 
detections northwards of the optical neb- 
ula, although there is neither an enhance- 
ment nor a deficit of stars at the bright 
photodissociation region (PDR). Although 
the reduced sensitivity at the chip gaps 
between the ACIS-I CCDs is not taken 
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into account, we doubt that the number 

of missing sources in the chip gaps is large 
enough to alter the spatial distribution of 
the stellar cluster. 

Projected stellar surface density (num- 
ber of stars arcmin"^) contour maps for the 
unobscured and obscured populations are 
shown in Figure 11, in logarithmic scale. 
The stellar density is smoothed using a 0.5' 
radius kernel. As expected, the highest 
concentration of stars is at the core of the 
cluster, where the massive O stars are lo- 
cated. The distribution of stars is rather 
spherically symmetric centered at the 03 
stars, except for a density enhancement 
centered at the tip of the S.E.C. (centered 
at RA=17:24:48.3, Dec=-34: 15:05). Fig- 
ure 12 further elucidates the structure of 
the stellar distribution by showing the ra- 
dial profile of the stellar surface density. 
The stellar distribution appears to con- 
sist of two components: a compact clus- 
ter within the central 1.5 pc with a log- 
normal distribution, and a cluster halo ex- 
tending to 8 pc with an approximately ex- 
ponential distribution. The number of X- 
ray stars in the central region is 333 while 
the extended region contains 446 stars. Be- 
cause these stars are poorly studied, we 
can only speculate that the stars in the 
central region represent the coeval dense 
cluster core of Pismis 24, while the appar- 
ent extended wing is a mixture of evapo- 
rating cluster members and newly formed 
stars distributed throughout the molecular 
cloud enveloping the cluster. 

Figure 12 also compares the radial pro- 
file to the profile of the Orion Nebula 
Cluster within ~1.5 pc for Chandra stars 
(Getman et al. 2005a; Feigelson et al. 
2005). The ONC is only extensively stud- 



ied within a projected radius of ~2.5 pc 
from 6^ Ori C (sec Hillenbrand & Hart- 
mann 1998 for a review of ONC struc- 
ture and dynamics). It is strongly concen- 
trated, with a central density > 200 times 
that of the widespread cloud stellar popu- 
lation distributed at radii > 2.5 pc. The 
ONC could possess an extended halo of 
stars as we see in the more massive Pismis 
24 cluster but it would be confused with 
the surrounding stellar population. Thus 
we limit the comparison of the two X-ray 
selected clusters to radii <1.5 pc. 

The stellar density of the ONC core ap- 
pears to be ~5 times higher than the X- 
ray population in the Pismis 24 cluster 
core, but this can be largely attributed 
to the differenccin the number of faint 
stars detected: the ONC sample from the 
800 ksec COUP observation extends into 
the brown dwarf regime, while the Pis- 
mis 24 sample becomes incomplete around 
M < 3M0 (§5.1). Extrapolating to brown 
dwarf masses using the COUP mass distri- 
bution (Getman et al. 2005a), we estimate 
that the NGC 6357 cluster central density 
would be ~20 times higher than the cur- 
rently observed one if it were as sensitive as 
the COUP observation. We thus conclude 
that the NGC 6357 cluster central density 
is probably ~4 times higher than the ONC 
central density and similar to the central 
density of the NGC 3603 super star cluster 
(sec footnote 1 in Hillenbrand & Hartmann 
1998). 

Figure 13 zooms in to the core region 
of the cluster using the SIRIUS K im- 
age as background. The ONIR identi- 
fied X-ray sources (circles) and unidenti- 
fied X-ray sources (crosses) are shown with 
the known massive members highlighted 
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in white. Only a few of the unidentified 

sources appear isolated, without infrared 
sources nearby. Most are clustered in the 
15" neighborhood of the two 03 stars; 24 
are located ~ 1.5—2" from near-lR sources, 
therefore the available NIR images may not 
be able to resolve them. These unidenti- 
fied sources very likely are real low-mass 
members concentrated around the massive 
stars. 

5.3. The Morphology of the Ring- 
like Nebula 

A puzzling feature arises when we com- 
pare the core of the massive cluster cen- 
tered on the two 03 stars to the large 
ring-like nebula with a large cavity that 
makes up the wider view of NGC 6357 
(Figure 1): the cluster is significantly off- 
center to the north. The "diamond ring" 
morphology of NGC 6357 is apparent in 
the Ha and mid-lR images, with the bright 
point along the edge brightened torus. In 
contrast, the Rosette Nebula has a sim- 
ilar annular nebula morphology, but the 
massive OB stars in NGC 2244 exciting 
the nebula are located right in the cen- 
ter. The hollow Hll morphology is at- 
tributed to strong OB stellar winds exca- 
vating molecular materials and the deple- 
tion of gas by newly formed stars (Dor- 
land et al. 1986; McCuUough 2000; Towns- 
ley et al. 2003). Therefore the nebular 
shape is generally accepted as spherical or 
cylindrical. In NGC 6357, the most mas- 
sive members are so close to the north- 
ern edge of the nebula that it seems im- 
plausible that the cavity is created by the 
strong winds of the current generation of 
hot stars. In general, supernovae explo- 
sions and strong stellar winds can shape 



the natal cloud by creating wind-blown 
"bubbles" -cavities. However, there is no 
evidence for a recent supernova: neither 
diffuse X-ray emission from a hot tenuous 
plasma nor non-thermal radio emission is 
detected in this region. We further note 
that the presence of WR 93, an evolved 
star, is circumstantial evidence that some 
population of massive stars in NGC 6357 
is old enough to produce supernovae. 

It is possible that a previous-generation 
of hot stars created the large cavity and 
formation of the current ionizing popula- 
tion was triggered at the densest edge of 
the excavated cavity. Combined with dif- 
ferent viewing angles, this can result in di- 
verse morphologies such as an egg, or a 
diamond ring, rather than a symmetrical 
annulus. Ha images of quite a number of 
regions where the stars are displaced with 
respect to the center of an evacuated cav- 
ity have been shown in the Large Magel- 
lanic Cloud (Gaustad et al. 1999) and mod- 
eled as modified Stromgren spheres (Mc- 
CuUough 2000). The key for testing this 
hypothesis is to identify the older cluster. 
X-ray sampling provides a critical test as 
X-ray luminosities decay only slightly over 
~ 10 Myr (Preibisch & Feigelson 2005). 
The sub cluster seen in the surface density 
map ~ 4' southeast of the Pis 24 cluster 
could be the remnant core of the previ- 
ous generation of stars. ONIR photometric 
and spectroscopy study of the ~30 stars in 
this subcluster, permitting their placement 
on the HR diagram, might reveal whether 
they are older than the main cluster con- 
centrated off the center of the nebula. 

An alternative scenario is that the cav- 
ity is indeed created by the massive Pis 24 
stars that are emitting ionizing photons. 
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In that case, the displacement from the 
center could be explained with spatial pro- 
jection and an inhomogeneous molecular 
cloud, as first proposed by Bohigas et al. 
(2004). The O stars were close to the geo- 
metrical center of the Pismis 24 cluster and 
created an expanding HII bubble, which 
may have encountered a much denser inter- 
stellar medium in the northern part than 
in other directions. The apparent deficit 
of X-ray sources to the north may be evi- 
dence for the absence of star formation in 
the densest region, where a new generation 
of stars may emerge eventually. 

6. X-rays across the Mass Spectrum 

6.1. X-rays from known massive 
stars 

OB stars have been known to be X-ray 
sources since early observations from the 
Einstein satellite (Harnden et al. 1979; 
Seward et al. 1979). Models were de- 
veloped where instabilities in radiatively- 
driven stellar winds from massive stars 
produce shocks and heat the gas to X- 
ray emitting temperatures (Lucy & White 
1980; Owocki et al. 1988; Owocki & Co- 
hen 1999). However, recent studies indi- 
cate that more complex models are needed 
to account for unexpected X-ray emission 
line profiles and hard, variable continuum 
emission (e.g., Waldron et al. 2004; Stelzer 
et al. 2005). 

We detect X-ray emission from all 
known early type stars ranging from 03 
to BO. 5 in this 38 ks ACIS observation. 
Table 6 summarizes ACIS detections and 
non-detections of X-ray emission from O 
and early B stars in the Pismis 24 cluster 
that are classified in literature. 



Two of the newly classified 03 stars 
Pis 24-1 (03.51f) and Pis 24-17 (03.51111; 
Walborn et al. 2002) exhibit X-ray lumi- 
nosities Lt^c ~ 10^^ ergs s^^ This is 
consistent with HD 93128 in the Carina 
Nebula cluster Tr 14, which has similar 
spectral type (Evans et al. 2003). Their 
lightcurves remain constant through the 
observation, which is typical for O stars, 
although variability in X-ray emission from 
O stars hke ( Ori is well known (Berghoe- 
fer & Schmitt 1994). Figure 14 presents 
their X-ray spectra and spectral fits (as- 
sumed abundances of 0.3 Zq unless other- 
wise noted). Compared to the best fit us- 
ing a single temperature plasma model, the 
X-ray spectrum of Pis 24-1 (^344) is bet- 
ter fit by a two temperature plasma model 
{Nh = 1.3 X 1022cm-2, /cTi = 0.5 keV, 
and kT2 = 1.7 keV). Pis 24-17 (#420) 
can be adequately fitted by a single soft 
component with Nh = 1.4 x lO^^cm^^ 
and kT = 0.7 keV. The fit is improved 
by adopting enhanced Mg (0.6 x Zq), Si 
(0.6 X Zq), and Ar abundances (2.7 x Zq). 
The Si XIII fine at 1.86 keV in Pis 24- 
17 is exceptionally strong, comparable to 
the same line seen in the O star IRS 2 
in RCW 38 (Wolk et al. 2002). The de- 
rived temperatures for soft and hard com- 
ponents for these O stars are similar to the 
X-ray temperatures of other single mas- 
sive stars: some O stars exhibit a simple 
single temperature plasma ~ 1 keV (e.g., 
MacFarlane et al. 1993; Rho et al. 2004; 
Townsley et al. 2006a), while a few re- 
quire an additional high energy component 
(e.g., Corcoran et al. 1994; Kitamoto & 
Mukai 1996; Evans et al. 2004; Stelzer et al. 
2005). Gagne et al. (2005) shows that the 
magnetically channeled wind shock model 
(Babel & Montmerle 1997a,b) with strong 
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line-driven winds in a single O star can 
adequately reproduce both the soft and 
the hard compents in the recent phase- 
resolved Chandra grating spectra of 9^ Ori 
C (05.5V). In some cases, the high energy 
component may in fact be an indication 
of binarity, as powerful winds in two mas- 
sive components shock to produce hard X- 
rays (e.g., Portegies Zwart et al. 2002; Al- 
bacete Colombo et al. 2003; Townsley et 
al. 2006b). 

The brightest X-ray point source in the 
ACIS field of view is the WC7+07-9 bi- 
nary WR 93 (#747) with a Lt,c = 1-6 x lO^^ 
ergs s~\ roughly consistent with the ear- 
her ROSAT PSPC value Lx = 1.25 x 10=^^ 
ergs s-^ (0.2-2.4 keV; Pollock et al. 1995). 
It is one of the spectroscopic WR bina- 
ries that are bright enough to be detected 
by previous generation X-ray telescopes; 
single WC stars are much fainter (Oski- 
nova et al. 2003). Thus the strong X-ray 
emission from binary WR stars is likely 
generated in the region where the mas- 
sive stellar winds collide. However, while 
archetype coUiding-wind binaries such as 
WR 140 (WC7+04V) display character- 
istic non-thermal emission in the radio 
band (Dougherty et al. 2005), WR 93 
shows thermal radio emission (Abbott et 
al. 1986). Also the X-ray hghtcurve ap- 
pears to be constant through the entire 
observation; no orbital X-ray variability is 
seen ^. The X-ray spectral fit for WR 93 
is shown in Figure 14 with an absorption 
column Nff — 7.4 x 10^^ cm~^ and a rather 
hard kT — 2A keV. The large absorption 
column densities towards the two 03 stars 
and the WR star derived from X-ray spec- 

^No orbital period for WR 93 is reported in litera- 
ture. 



tral fits match well with their large visual 
extinction {Ay ~ 5 — 7) and the reddening 
from JHK colors. 

Five of the bright Pismis 24 stars in 
Moffat & Vogt (1973) - Pismis 24-4, 7, 
8, 9, 11 - were not included in the spec- 
tral classification observation of Massey et 
al. (2001). Based on their strong X-ray 
emission (except for Pis 24-9 which is un- 
detected) , absorption column derived from 
X-ray spectral fitting, UBV brightness, 
and JHK colors, we suggest that they are 
also young OB stars and cluster members 
(see Table 7). 

X-ray emission from OB stars was re- 
ported from ROSAT observations to dis- 
play a characteristic efficiency Lx/Lioi ~ 
10~^ (Chlebowski et al. 1989; Berghoefer et 
al. 1997), and this empirical ratio has been 
examined in recent X-ray observations of 
massive stars (e.g., Corcoran 1999; Evans 
et al. 2004; Stelzer et al. 2005). Prom 
the COUP observation, the emission ratio 
for late O and B stars is found to scat- 
ter around the canonical Lx / Li,oi value by 
three orders of magnitude. However, as 
they cautioned (Stelzer et al. 2005), this 
may not represent the most luminous early 
O stars due to the absence of spectral types 
earlier than 07. Therefore the Pismis 24 
cluster, with two 03 stars, is particularly 
valuable for assessing this relation in a sin- 
gle massive cluster. Using Mboi reported in 
Massey et al. (2001) Figure 15 shows 

'It should be noted that Heap et al. (2006) have 

proposed a new calibration between the spectral 
type and effective temperature of O stars, based 
on analysis of high resolution HST, FUSE, and 
optical spectra. The spectral calibration in gen- 
eral use is by Vacca et al. (1996), which Massey 
et al. (2001) used to derive bolometric luminosity 
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the Lx-Liioi relation for our detections of 
cataloged O and BO stars together with 
the strong wind sample and weak wind 
sample from the COUP study (Stelzer et 
al. 2005). The Lx/L^oi ratio for O stars 
in the Pismis 24 cluster is consistent with 
the ~ 10"''' value with an order of mag- 
nitude scatter. Although there are only 9 
data points and one upper limit (L^o; for 
WR 93 is not well determined due to the 
poorly characterized bolometric correction 
for late-type WC stars; Massey et al. 2001), 
it appears that the Lx / L^oi correlation has 
less scatter for stars of the earliest spectral 
types and the physics behind this relation 
is hkely to be physically related to stellar 
wind properties. 

6.2. Newly-discovered Candidate O 
Stars 

We assembled a list of candidate O stars 
in Table 7 to facihtate future spectroscopic 
follow-up via two approaches. Table 7a 
lists the 13 sources brighter than 80 counts 
(to guarantee a reliable X-ray spectral fit) 
with 0.5-8.0 keV absorption corrected X- 
ray luminosities logL^^c > 32.0 ergs as 
potential new high-mass stars, since pre- 
MS stars rarely reach X-ray luminosities 
> 10=^2 g^gs g-i (Favata & Micela 2003; 
Getman et al. 2005b). Then we examined 

for Pismis 24 O stars. The derived log Lboi us- 
ing the Heap et al. (2006) cahbration is 0.4 dex 
lower than log Lboi using the Vacca et al. (1996) 
value for an 04V star, and 0.15 dex lower for an 
08V star. Here we use the new bolometric correc- 
tions from Heap et al. (2006) to derive bolometric 
luminosities for Pismis 24 stars. The bolometric 
correction value for the two 03.5 stars is not pro- 
vided by Heap et al. (2006), and we adopt —0.6. 
COUP data points are from Stelzer et al. (2005) 
and originally derived from Hillenbrand (1997). 



their X-ray light curves and observed opti- 
cal and infrared properties to exclude pos- 
sible AGNs. Since certain O stars are ob- 
served to exhibit powerful flares (Feigelson 
et al. 2002), we do not exclude the X-ray 
luminous sources that may be extraordi- 
narily bright flaring pre-MS stars. Their 
lightcurve variability is noted in Table 7. 

O stars may reveal themselves as lumi- 
nous in infrared but appear relatively faint 
in X-rays due to heavy obscuration in their 
early evolutionary stages. For example, an 
09V star emitting logL^ c = 32.0 ergs s~^ 
behind 15 mag of visual extinction would 
have K = 10 and ~ 40 counts in our ACIS 
image when observed at a distance of 2.5 
kpc (assuming a kT ~ 0.6 keV and with- 
out considering a possible K-cxccss). We 
selected bright near-lR sources (Kg < 10.0 
mag) detected in the X-ray image to ex- 
amine whether their JHK colors and X-ray 
properties are suggestive of obscured early 
type stars. This K brightness selection is 
used in Hanson, Howarth, & Conti (1997) 
when they select the spectroscopy sample 
of massive YSOs in M 17, and we further 
adapt it to X-ray emitting samples. All 
9 known members with spectral type ear- 
lier than BO. 5, the WR-l-07 binary, and 
the foreground A81V star were recovered 
using this critieria. Table 7b gives 11 ad- 
ditional candidates that have colors con- 
sistent with early-type stars and are not 
known foreground stars. Figure 7 and 
Figure 8 convincingly suggest that highly 
reddened bright infrared source #654 and 
#694 are indeed very obscured early O 
stars with some IR excess. Source #140 
is almost certainly a high-mass star since 
it qualifles using both the IR and the X-ray 
selection criteria. As a cautionary note for 
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this method, the presence of X-band ex- 
cess may mislead us to include some accret- 
ing intermediate-mass stars as candidate O 
stars. 

Given the uncertainties in distance, age, 
and absorption, readers are cautioned that 
ONIR spectroscopy on these objects (Ta- 
ble 7) should yield more appropriate classi- 
fications for them (e.g., Massey et al. 2001; 
Walborn et al. 2002). The X-ray sample 
suggests that the optical sample of Massey 
et al. (2001) may have significantly under- 
estimated the O star population. However, 
using a smaller distance of 1.7 kpc will re- 
sult in a 0.35 dex decrease in the derived X- 
ray luminosity. Some of the X-ray selected 
candidates may indeed be lower mass fiar- 
ing pre-MS stars, which are readily distin- 
guishable from their ONIR spectra. If a 
large fraction of these 24 candidate O stars 
are confirmed, the OB population of Pis 24 
(Table 6) is doubled to tripled. 

On the top right corner of the infrared 
color-magnitude diagram (Figure 7 and 
Figure 8), there is another mysterious 
bright source, #19 (CXOU J172413.60- 
341456.7), whose color and brightness 
{J - H = 2.006, H - K = 0.991, K = 
5.334) suggest that it is a very luminous 
object subjected to significant absorption 
{Ay ~ 15 — 20). It has 17 net counts 
in X-ray and the median photon energy 
is 2.5 keV, which is consistent with being 
obscured and not a foreground star. How- 
ever, if the star is located at 2.5 kpc, the 
intrinsic K-band magnitude and infrared 
luminosity will be exceptionally luminous. 
Thus it is plausible that source #19 is 
a massive YSO or a young post main- 
sequence star like WR 93. Spectroscopy 
is highly warranted and should reveal the 



nature of this unusually luminous but ob- 
scured star. 

6.3. Intermediate Mass Stars 

The large population of the Pismis 24 
cluster offers an exceptionally rich sam- 
ple of intermediate mass stars. Of par- 
ticular interest are those with infrared ex- 
cess, possibly due to the presence of pro- 
toplanetary disks. We note that the es- 
timated stellar masses from the infrared 
color-magnitude diagrams rely on the un- 
certain age and distance, which are de- 
termined with improved precision in re- 
cent literature. The ambiguity is more 
severe among the early B stars (Figure 7; 
for example, the inferred stellar mass for 
a 14 Mq star from the 1 Myr isochrone 
can be as low as 4 Mq when estimated 
from a 0.3 Myr isochrone. However, the 
X-ray selected sample is important for fu- 
ture follow-up for this poorly studied re- 
gion. We hst in Table 8 100 X-ray selected 
candidate intermediate-mass stars, which 
have dcrcddencd spectral types AO— BO 
estimated from the NIR color-magnitude 
diagram (Figure 7 and Figure 8) assum- 
ing 1 Myr for age and 2.56 kpc for dis- 
tance. Optical spectroscopy is needed 
to give accurate spectral classifications. 
Among them, four stars exhibit significant 
X-band excesses in the NIR color-color 
diagram (Figure 6), which are good can- 
didates for HAeBe stars. Ha spectroscopy 
on this sample might further reveal the ex- 
istence of accretion. Their X-ray luminosi- 
ties are on the order of log Lx ~ 30.0—31.5 
ergs s~^, consistent with previous ROSAT 
PSPC (0.1-2.4 keV) observation of HAeBe 
stars (Zinnecker & Preibisch 1994). 

One interesting result has been men- 
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tioned in §4.5 - the available IR photome- 
try data suggest that the ratio between the 
number of Class 11 and Class 111 pre-MS 
stars with mass > 2Mq is low, ~ 1 : 15. 
Over the mass range 2 — 16Mq, the ra- 
tio between the number of Class II and 
Class 111 pre-MS stars is even lower, ~ 4%. 
It is important to note that such cl Tel- 
tio is derived from X-ray selected stars, 
without any prior knowledge of disk in- 
dicators. The low fraction of the cluster 
members possessing disks is rather surpris- 
ing since apparently these X-ray emitting 
pre-MS stars are very young (in the com- 
pany of 03 stars). However, the low disk 
frequency is consistent with the previous 
studies in NGC 6611 by Hillenbrand et al. 
(1993). They find that optically thick cir- 
cumstellar disks are already rare among 
the intermediate-mass pre-MS stars with 
ages less than 1 Myr and suggest the disk 
lifetimes are much shorter for the massive 
stars than those of solar type stars. It 
seems also consistent with the suggestion 
that disk evolution happens rapidly in clus- 
ters based on recent Spitzer observations 
(e.g., Hartmann et al. 2005; Sicilia-Aguilar 
et al. 2006; Lada et al. 2006). The dras- 
tic radiation environment may play a big- 
ger role here in the dissipation of disks, 
rather than mass loss via accretion. Fur- 
ther Spitzer study on this statistically sig- 
nificant sample will provide constraints on 
the timescale of the circumstellar disk dis- 
sipation and the formation of planetary 
systems. 

6.4. Flaring pre-MS Stars 

There are 31 sources identified as highly 
variable {Pks ^ 0.005; variability fiag "c" 
in Table 1 and Table 2); their lightcurves 



exhibit flaring activity as most frequently 
seen in pre-MS stars when magnetic re- 
connection events occur in their coronae. 
Eleven of them have more than 80 net 
counts (0.5-8.0 keV), so we can rehably de- 
rive their X-ray luminosities through spec- 
tral fits. Seven of the fiaring sources ex- 
hibit luminosities above logL^^c > 32.0 ergs 
s~^, which is exceptionally high for pre- 
MS stars in their flaring phase (Feigelson 
& Montmerle 1999; Favata & Micela 2003). 

One extraordinarily intense fiare is 
seen in source #672 (CXOU J172457.87- 
341203.9), comparable to some of the 
strongest X-ray fiares known in pre-MS 
stars (Grosso et al. 2004; Favata et al. 2005; 
Getman et al. 2006). The time-energy di- 
agram, binned lightcurve, and X-ray spec- 
trum for #672 are shown in Figure 16. In 
contrast to the typical impulsive fast rise 
phase of 2 hr observed in the X-ray flares of 
pre-MS stars, this flare shows a ~ 4 hr ris- 
ing phase to reach its highest fiux and still 
remains at the peak level at the end of our 
observation. This is similar to the pow- 
erful flare seen in the non- accreting pre- 
MS star LkHcK 312 in the Orion B cloud 
(Grosso et al. 2004). The count rate is 
~ 15 times higher during flaring compared 
to the quiescent level before flaring. The 
spectrum is hard {{E) — 3.0 keV), as seen 
in many pre-MS stars in COUP during 
flares (Favata et al. 2005). The spectral 
fitting gives a rather hard kT > 10 keV, 
corresponding to thermal plasma of > 120 
MK. The infrared countpart for this X-ray 
source shows K — 10.2 and its location 
in the color- magnitude diagram suggests a 
high-mass star. The SIRIUS images were 
visually examined and there is a possibil- 
ity that the flare originates from a poorly 
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resolved low mass companion to the high- 
mass star. If a future spectrum confirms 
that it is an early type star, the flaring 
behavior is similar to that observed in the 
Orion Nebula 09.5 star 0^ Ori A (Feigel- 
son et al. 2002), either from an extraordi- 
nary magnetic reconnection flare in a low 
mass companion, or from non-standard 
wind shocks in the massive stellar winds 
(Feigelson et al. 2002; Stelzer et al. 2005). 

6.5. X-ray Selected Deeply Embed- 
ded Population 

In Table 9 we assemble a sample of 
16 heavily obscured objects. These are 
sources that have (E) > 3.0 keV (with > 
20 net counts to be rehable) or log Nh,x > 
22.5 cm~^ derived from spectral fitting. 
This subsample of the hardest sources 
(classified H in Table 5) likely contains 
deeply embedded protostars in addition to 
older stars heavily obscured by molecu- 
lar material. Their spatial distribution is 
shown in Figure 17, overlaid on a DSS op- 
tical image. Several of them appear to be 
clustered in the optical/infrared dark col- 
umn (S.E.C. in the CO map) clearly seen 
in visible and 2MASS images (Figure 1). 
Others are widely distributed along the 
rim of the ringlike nebula. A plausible ex- 
planation is that the nebulosity traces the 
rim of a bowl-shape shell facing us, thus 
the column density in our line-of-sight to 
those objects is much higher. This geo- 
metric configuration is consistent with the 
previous conclusion that G353.2-I-0.9 is a 
bhster HII region viewed face-on (Massi et 
al. 1997). 



7. Diffuse X-ray Emission 

With two of the most massive stars 
known and a luminous WR star in the field, 
it is reasonable to expect diffuse X-ray 
emission in NGC 6357, from strong stel- 
lar winds and parsec-scale wind-wind colli- 
sions (Townsley et al. 2003). A more recent 
study of X-ray diffuse emission in the 30 
Doradus star-forming complex is reported 
in Townsley et al. (2006b), which provides 
our technical procedures to extract diffuse 
emission. From the point source locations 
and extraction regions (Table 1 and 2), 
we can mask the point sources and ex- 
tract the spectrum of possible remaining 
diffuse emission. Two regions of different 
scales are used for extraction, which yield 
consistent measurements of the emission. 
They are outhned by the white and the 
cyan polygons in Figure 2b. The back- 
ground spectrum was scaled by the effec- 
tive area and by the ratio of the geometric 
areas to be appropriate for the region of 
diffuse emission being fit; see Section 5.1 
of Townsley et al. (2003) for details. The 
spectra can be fit using a soft optically 
thin thermal plasma with kT ~ 0.6 keV 
and logiVi^ ~ 22.1 cm~^. The correspond- 
ing observed soft band X-ray luminosity is 
logL^ ~ 32.4 ergs s~^ (logL^ c ~ 33.6 ergs 
s~^ if corrected for absorption). This is 
comparable to the observed diffuse X-ray 
emission in M 17 and the Rosette Nebula 
(Townsley et al. 2003). 

However, the detected emission does 
not necessarily arise from wind-generated 
plasma. Under the assumption that the 
gas temperature within the cluster is rela- 
tively constant and hot (> 10^ K), Stevens 
& Hartwell (2003) show that the luminos- 
ity level of the expected diffuse emission in 
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a stellar cluster can be estimated using a 
scaling factor (equation 10), which Getman 
et al. (2006) adopt to evaluate the nature 
of X-ray diffuse emission in the Cepheus 
B region. We follow the same method to 
obtain a crude estimate of diffuse emission 
from the massive star winds Two 03 
stars dominates the cluster core and two 
parameters are needed to predict the X- 
ray luminosity: the total stellar mass loss 
rate M*, and the mean weighted termi- 
nal velocity for the stars in the cluster v^,. 
For the 03 star Pis 24-1, Benaglia et al. 
(2001) derived a terminal wind velocity of 
~ 2295 km s~^ and an upper limit for 
the mass loss rate of M < 4.5 x lO^^M© 
yr~^. No measurements for Pis 24-17 is 
found in the literature, and we assume 
Voo and M are comparable to Pismis 24-1. 
Adopting these wind parameters, and an 
unreahstically optimistic 100% efficiency 
for conversion from stellar wind to ther- 
malization of the cluster, and ~ 2 pc for 
the radius of the cluster core, we estimate 
the diffuse X-ray luminosity Lx < 10^^ 
ergs s~^ (equation 8 and 10). This is ap- 
parently insufficient to account for the ob- 
served luminosity. 

An alternative explanation is that the 
diffuse emission is mainly the contribution 
from unresolved X-ray emission from thou- 
sands of low mass pre-MS members in Pis- 
mis 24. Assuming that Pismis 24 and the 
ONC have identical IMFs with identical 
XLFs, the only difference here is that Pis- 
mis 24 cluster has 5 times the total known 
population of the ONC (see §5.1). We 

The WC+0 star WR 93 is not in the cluster core 
region, thus we do not include it in the calculation. 

For reader's reference, the mass loss rate and the 
terminal wind velocity for WR 93 is given in §2 



focus on the lightly obscured populations 
since they are less affected by the different 
amount of absorption along the two sight- 
lines. Scaling from COUP (Feigelson et al. 
2005), the total expected soft-band X-ray 
luminosity is log L^^c ~ 33.8 ergs s~^ from 
the entire Pismis 24 population and the ac- 
tual observed soft luminosity from resolved 
cluster members is logL^ ~ 33.2 ergs s~^ 
(since absorption correction to individual 
stars is infeasible, we consider this to be 
a lower limit to the intrinsic soft-band lu- 
minosity Lg^c)- Thus the integrated X-ray 
emission from the unresolved X-ray emit- 
ting pre-MS population in Pismis 24 is es- 
timated to be log Ls^c < 33.7 ergs s~^, af- 
ter deducting the observed Ls from the ex- 
pected Ls.c- This is fully consistent with 
the absorption-corrected soft diffuse X-ray 
luminosity log L^ c ~ 33.6 ergs s~^. There- 
fore we suggest that most of the diffuse 
emission is probably the combined con- 
tribution from individually undetected X- 
ray-faint pre-MS cluster members. 

We further note that X-ray emission is 
clearly present on the ACIS-S chips. Sev- 
eral point sources are detected, with very 
large PSFs due to large off-axis angles. 
The X-ray emission captured at the lower 
edge of the S-array field corresponds to the 
location of the H ll region G353.2-I-0.7 (see 
Figure 2b and Figure If). But due to the 
poor PSF, no further conclusion can be 
drawn on whether the emission is diffuse 
in nature. Nevertheless, this detection con- 
ffims the ROSAT detection of G353.2+0.7 
and implies that it contains lots of X-ray- 
emitting young stars. 
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8. Summciry 

This 38 ks ACIS observation of tlie NGC 
6357 field fias provided tfie first fiigli spa- 
tial resolution X-ray image of this massive 
star-forming complex. We summarize the 
main results of our study as follows: 

1. We detect 779 X-ray sources with 
a limiting X-ray sensitivity of ~ 10^° ergs 
s~^. There are 445 ACIS sources matched 
to 2MASS point sources with < 2" sep- 
aration. Further visual inspection of an 
archival HST image and SIRIUS JHK im- 
ages yield a total of 616 ONIR counter- 
parts. We estimate the extragalactic con- 
taminants via a careful simulation of de- 
tection efficiency of the background AGN 
population, and conchidc that no more 
than 4% of our detections are AGNs. Sim- 
ilarly we adopt the Besangon stellar popu- 
lation sythesis model of the Galactic disk 
and suggest 11 stars as foreground ob- 
jects based on their X-ray softness and 
NIR color. Excluding ~ 20 X-ray sources 
without ONIR identifications that are pos- 
sible AGNs, the rest of the ~ 140 X- 
ray sources without ONIR counterparts are 
likely new Pismis 24 members that arc 
clustered around the massive stars and/or 
deeply embedded in the cloud. The X-ray 
detected population provides the first deep 
probe of the rich population of this massive 
cluster, increasing the number of known 
members from optical study by a factor of 
~ 50. 

2. On the NIR color-color diagram (Fig- 
ure 6), it is clear that ACIS detects more 
Class 111 sources than Class II/I sources. 
In the color-magnitude plot (Figure 7), 
the locations of most of the ACIS sources 
are consistent with being high mass stars 



and pre-MS intermediate mass stars (M > 
2M0), subjected to a visual extinction of 
5 < Ay < 15 at 2.5 kpc. The new high 
mass stars, if spectroscopically confirmed, 
double the known high mass population of 
the cluster. Dozens of intermediate mass 
stars are newly recognized members of the 
cluster. The ratio between Class III and 
Class II objects in the intermediate- to 
high-mass range is ~ 15 : 1. This is 
an important new sample for stellar and 
disk studies of intermediate mass pre-MS 
stars. The brightness of several K-excess 
sources makes them extremely favorable 
for ground based follow-up. We find a very 
luminous X-ray emitting infrared source 
subjected to significant absorption, which 
may be a new obscured massive YSO or a 
Wolf-Rayet star. 

3. We determine the spatial distribution 
of over 700 new members to sub-arcsecond 
precision for the first time. The cluster 
is roughly spherically symmetric, with the 
highest stellar density around the massive 
03 stars. A density enhancement centered 
at the tip of the molecular filament de- 
serves further investigation (Figure 11). A 
radial profile of the stellar surface density 
shows that the physical size of the cluster 
spans over 2 pc in radius and extends be- 
yond 6 pc (Figure 12). 

4. XLFs are derived for lightly obscured 
and heavily obscured stellar populations in 
NGC 6357 (Figure 9). The high X-ray lu- 
minosity end (logLx > 30.3 ergs s^^) of 
the NGC 6357 XLF is clearly consistent 
with a power law relation as seen in COUP 
and Cep B, allowing a scaling comparison 
with the COUP XLF. The first estimate of 
the total cluster population of Pismis 24 is 
a few times the known Orion population. 
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The presence of two 03 stars is consistent 
with the standard IMF. 

5. We detect all 10 known OB stars 
with spectral type earlier than Bl in the 
Pismis 24 cluster, including the WC7+07 
binary WR 93. Their Lx/Ljjoi ratios are 
consistent with the canonical 10~^ value 
with less than an order of magnitude scat- 
ter. The 03 stars Pis 24-1 (03.5If) and 
Pis 24-17 (03.5IIIf) exhibit X-ray lumi- 
nosities Lx ~ 10^^ ergs s~^. The X-ray 
spectrum of Pis 24-1 shows a soft emission 
component {kT ~ 0.5 keV) and a harder 
component {kT ~ 1.7 kcV), while Pis 24- 
17 exhibits an unusaily strong Sixill line 
with a soft 0.7 keV thermal plasma. The 
WC7-F07 binary WR93 is the brightest X- 
ray source in the field (Lx — 2 x lO'^^). 
Its hard spectrum suggests that the X-ray 
emission is generated from colliding winds. 
The hghtcurves of these massive stars are 
constant through the observation. We as- 
sembled a list of candidate O stars from 
sources with high X-ray and infrared lu- 
minosity to facilitate future spectroscopy 
follow-up. Four candidate HAeBe stars are 
also identified. 

6. We report the detection of X-ray 
emission from an EGG at the tip of a 
molecular pillar, which was previously 
found to be an IR source and the peak 
of the radio continuum. It is in the early 
evolutionary of stages of star formation. 
The spectral type estimated from its SED 
is B0-B2. The non-detection of EGGs in 
the Eagle Nebula may be the result of 
lower stellar mass and earlier evolution- 
ary phases. We assemble a sample of 16 
deeply embedded objects. 

7. Several flaring sources exhibit lu- 
minosities above log L^^c > 32.0 ergs s~^, 



which is exceptionally high for pre- MS 
stars. One powerful flare seen in ACIS 
#672 (CXOU J172457.87-341203.9) is 
comparable to the strongest X-ray flares 
known in pre-MS stars (L^c ~ 32.5 ergs 
s"-*^). Its lightcurve rises 15 times above 
the quiescent level during the flare, with 
a slow rise phase that is distinct from the 
fast rise flares seen pre-MS stars. 

8. Soft unresolved X-ray emission in the 
NGC 6357 region is present. However, the 
luminosity derived from spectral fitting is 
consistent with the estimated level of in- 
tegrated emission from the unresolved pre- 
MS stars. 
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Table 1 

Main Chandra Catalog: Basic Source Properties 



Source Position Extracted Counts^ Characteristics 



Seq 


CXOU J 




*J2000 


Err 


9" 


Not 


ANet 


Bkgd 


Not 


PSF 


Signif 


log Pb" 


Anom° 


Var"-' 


EffExp 


Mod B 


# 




(deg) 


(deg) 


(") 


{') 


Full 


Full 


Full 


Hard 


Frac 










(ks) 


(koV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


1 


172353.98-340850.8 


260.974932 


-34.147448 


0.4 


10.7 


267.2 


17.7 


22.8 


160.0 


0.91 


14.7 


<-5 




a 


32.7 


2.4 


2 


172358.42-340802.6 


260.993436 


-34.134056 


0.5 


10.2 


107.0 


11.5 


12.0 


71.7 


0.91 


8.9 


<-5 




b 


33.1 


2.6 


3 


172369.11-341217.6 


260.996328 


-34.204895 


0.8 


9.1 


23.2 


6.6 


12.8 


12.4 


0.91 


3.2 


<-5 




a 


33.0 


2.0 


4 


172401.20-340928.0 


261.005039 


-34.157783 


0.9 


9.1 


13.7 


5.3 


8.3 


7.2 


0.89 


2.4 


-4.3 




a 


33.7 


2.1 


5 


172402.27-341110.4 


261.009476 


-34.186237 


0.8 


8.5 


19.4 


5.8 


7.6 


20.0 


0.90 


3.1 


<-5 




a 


34.2 


3.7 


6 


172403.19-341402.2 


261.013313 


-34.233962 


1.1 


8.5 


16.6 


5.7 


9.4 


8.0 


0.90 


2.7 


<-5 




a 


34.2 


1.9 


8 


172405.42-341257.5 


261.022606 


-34.215981 


0.4 


7.9 


63.2 


8.9 


6.8 


58.3 


0.89 


6.7 


<-5 




a 


34.7 


3.9 


9 


172405.64-340710.0 


261.023532 


-34.119459 


0.7 


9.3 


38.2 


7.3 


7.8 


1.4 


0.90 


4.8 


<-5 




a 


33.6 


1.1 


10 


172406.70-341305.0 


261.027933 


-34.218067 


0.8 


7.6 


11.0 


4.8 


7.0 


9.3 


0.90 


2.1 


-3.5 




a 


34.9 


3.8 


12 


172407.65-341751.9 


261.031884 


-34.297766 


0.8 


9.3 


27.8 


7.1 


13.2 


13.7 


0.90 


3.7 


<-5 




a 


31.6 


2.1 



^Full band = 0.5-8 keV; Hard band = 2-8 keV. 
^Off-axis angle. 

Pb — probability extracted counts (full band) are solely from background. 

"^Source anomalies: g — fractional time that source was on a detector (FRACEXPO from mkarf) is < 0.9 ; c — source on field edge; p — source piled up; s — source on readout streak. 
Note that a reduced PSF fraction {significantly below 90%) indicates that the source is in a crowded region. 

^Source variability: a = no evidence for variability; b = possibly variable; c = definitely variable. No test is performed for sources with fewer than 4 total full-band counts. No value is 
reported for sources in chip gaps or on field edges. 

Note. — See §3.3 for a detailed description of the columns. The full table of 665 Chandra sources is available in the electronic edition of the Journal. 



Table 2 
Tentative Source Properties 



Source Position Extracted Counts^ Characteristics 



Seq 


CXOU J 


«J2000 


^J2000 


Err 




Not 


ANot 


Bkgd 


Net 


PSF 


Signif 


log Pb" 


Anom'^ 


Var° 


EffExp 


Mod E 


# 




(deg) 


(deg) 


(") 


(') 


Full 


Full 


Full 


Hard 


Frac 










(ks) 


(kcV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


m 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


7 


172403.64-340634.7 


261.015178 


-34.109657 


1.4 


10.0 


10.4 


5.1 


9.6 


2.6 


0.91 


1.8 


-2.7 




a 


32.4 


1.6 


11 


172406.84-341149.7 


261.028513 


-34.197157 


1.2 


7.5 


6.4 


3.7 


3.6 


4.0 


0.90 


1.5 


-2.4 




a 


33.0 


2.0 


16 


172412.29-341600.3 


261.051238 


-34.250107 


1.0 


7.0 


6.8 


3.9 


4.2 


3.7 


0.90 


1.5 


-2.4 




a 


33.1 


2.1 


23 


172415.22-340730.0 


261.063419 


-34.125011 


0.9 


7.5 


6.2 


3.6 


3.8 


3.8 


0.90 


1.3 


-1.8 




a 


32.2 


2.6 


24 


172415.43-340958.2 


261.064294 


-34.166176 


0.9 


6.2 


3.7 


3.0 


2.3 


0.5 


0.90 


1.0 


-1.5 






33.6 


1.4 


25 


172415.57-342044.7 


261.064909 


-34.345752 


1.5 


10.3 


7.9 


5.3 


13.1 


1.2 


0.91 


1.4 


-1.6 




a 


28.1 


1.4 


27 


172415.91-341441.2 


261.066326 


-34.244799 


0.8 


6.2 


6.3 


3.5 


2.7 


6.2 


0.90 


1.5 


-2.7 




a 


34.9 


3.5 


32 


172418.05-340824.5 


261.075226 


-34.140162 


0.9 


6.5 


4.6 


3.2 


2.5 


5.9 


0.89 


1.2 


-1.8 




a 


35.3 


4.5 


33 


172418.06-341650.6 


261.075253 


-34.280743 


1.0 


7.0 


6.6 


3.7 


3.5 


2.2 


0.89 


1.5 


-2.5 




a 


34.7 


1.1 


37 


172419.79-341130.3 


261.082481 


-34.191754 


1.0 


4.9 


4.8 


3.0 


1.2 


2.4 


0.90 


1.3 


-2.8 




a 


36.4 


1.1 



^Full band = 0.5-8 keV; Hard band = 2-8 keV. 
'^Off-axis angle. 

Pb ~ probability extracted counts (full band) are solely from background. 

"^Source anomalies: g — fractional time that source was on a detector (FRACEXPC) from rnkarf) is < 0.9 ; c — source on field edge; p — source piled up; s — source on readout streak. 
Note that a reduced PSF fraction (significantly below 90%) indicates that the source is in a crowded region. 

^Source variability: a = no evidence for variability; b = possibly variable; c = definitely variable. No test is performed for sources with fewer than 4 total full-band counts. No value is 
reported for sources in chip gaps or on field edges. 

Note. — See §3.3 for a detailed description of the columns. The full table of 114 Chandra sources is available in the electronic edition of the Journal. 



Table 3 

X-RAY Spectroscopy for Brighter. Sources: Thermal Plasma Fits 





Source 




Spectral Fit* 






X-ray Luminosities^ 




Goodness 


Notes 


Seq 


CXOU J 


logAfff 


kT 


log EM 


log is 


logi), logi),j, logLt 


log Lt,o 


of Fit 




# 
(1) 


(2) 


(cm-2) 
(3) 


(keV) 
(4) 


(cm-3) 
(5) 


(6) 


(ergs s~^) 
(7) (8) (9) 


(10) 


(11) 


(12) 



1 


172353. 


.98-340850. 


.8 


-0.08 


22. 


.0 


-HO. 


.10 




> 


10. 







-0.05 


55, 


.0 


-HO, 


.05 


31, 


.16 


31. 


.97 


32. 


.00 


32, 


.03 


32 


.18 


0, 


.92 


S; AGN? 


2 


172358. 


.42-340802. 


.6 


-0.2 


22. 


.3 


+0. 


.2 


-1, 


,1 


3.0 


-H2, 


.8 


-0.2 


54, 


.9 


-HO, 


.3 


30, 


.68 


31. 


.54 


31, 


.65 


31, 


.60 


31. 


.97 


0, 


.93 


2MASS 


51 


172423 


.93-340816. 


2 


-0.3 


21 


.4 


+0. 


.3 


-3, 


,7 


7.1 


-H2 


.9 


-0.09 


54. 


.3 


+0. 


,13 


30. 


.76 


31 


.22 


31 


.23 


31 


.35 


31 


.44 


1 


.07 


USNO 0558-0532169-H2MASS 


77 


172428 


.94-341450. 


.6 


-0.10 


21 


.5 


+0. 


.09 


-0, 


,8 


4.1 


-HI 


.2 


-0.05 


54. 


.9 


+0. 


,05 


31 


.40 


31 


.76 


31 


.77 


31 


.91 


32. 


.04 


1 


.06 


Pis24-15 


78 


172429 


.03-341813. 


.5 


-0.2 


21 


.9 


+0. 


,2 


-0, 


,8 


2.7 


-HI. 


.3 


-0.2 


54. 


.7 


+0. 


,2 


30. 


.88 


31. 


.32 


31. 


.37 


31 


.46 


31 


.72 


0, 


.27 


USNO 0556-0531216-H2MASS 


112 


172433. 


.04-341654. 


.5 


-0.1 


22. 


.1 


+0. 


.1 


-1, 


,1 


3.6 


-H2. 


.0 


-0.1 


54. 


,8 


-HO. 


,2 


30. 


.84 


31. 


.54 


31. 


.60 


31, 


.62 


31. 


.88 


0, 


.46 


USNO 0557-0528440-H2MASS 


140 


172434. 


.79-341318. 


.0 


-0.09 


21 


.9 


+0. 


,08 


-1 


,9 


7.0 


-H3. 


.0 


-0.05 


55. 


,0 


-HO. 


,06 


31. 


.27 


31. 


.96 


32. 


.00 


32. 


.04 


32. 


.20 


0, 


.71 


Pis24-ll 


158 


172436 


.63-341550. 


7 


-0.1 


22 


.0 


+0. 


.1 


-0 


,8 


2.9 


+ 1 


.2 


-0.1 


55 


,0 


+0. 


,2 


31 


.05 


31 


.63 


31 


.69 


31 


.73 


32 


01 


0. 


.80 


USNO ()557-0528485-H2MASS 


188 


172438 


.73-341202. 


.9 


-0.1 


22 


.0 


+0. 


.1 


-0 


,9 


3.:5 


+ 1 


,3 


-0.1 


54. 


8 


+0. 


,1 


30. 


9fi 


31 


.54 


31 


60 


31 


64 


31 


89 





66 


USNO ()557-0528507-H2MASS 


196 


172438 


.98-341220. 


6 


-0.1 


22 


.0 


+0. 


.1 


-0 


,B 


2.2 


+ 1 


,0 


-0.2 


54. 


,7 


+0. 


,2 


30. 


76 


31 


.21 


31 


.28 


31 


.34 


31 


.68 


0. 


.47 


USNO 0557-0528516-H2MASS 


203 


172439 


.51-341317. 


.9 


-0.2 


21 


.8 


+0. 


.2 


-1 


,7 


4.8 


+4 


,9 


-0.1 


54. 


,5 


+0. 


,1 


30. 


.79 


31 


.34 


31 


.37 


31. 


.45 


31. 


.61 


0. 


.69 


USNO 0557-0528522-H2MASS 


227 


172440 


.67-341403. 


.8 


-0.08 


22 


.1 


+0. 


.08 


-2 


,8 


7.6 


+2 


,4 


-0.07 


54. 


,9 


+0. 


,09 


30. 


.90 


31 


.81 


31 


.87 


31. 


.86 


32. 


.06 


0. 


.76 


2MASS 


243 


172441. 


.23-341138. 


.6 


-0.2 


22 


.0 


-1-0. 


.2 


-1, 


,0 


3.3 


-H2. 


,3 


-0.2 


54. 


,8 


-HO. 


,2 


30. 


.91 


31. 


.50 


31. 


.55 


31. 


.60 


31. 


.85 


1. 


.39 


USNO 0558-0532303-H2MASS 


270 


172441. 


.92-341158. 


.0 


-0.1 


22. 


.0 


-HO. 


.1 


-0, 


,6 


2.2 


-HO. 


,6 


-0.1 


66. 


,0 


-HO. 


,2 


31. 


.06 


31 


.50 


31. 


.66 


31. 


.63 


31. 


.96 


1. 


.31 


USNO 0558-0532308-H2MASS 


284 


172442. 


.14-341155. 


.9 


-0.2 


22. 


.0 


-HO. 


.2 


-3, 


,2 


6.9 


-H3. 


,1 


-0.10 


64. 


,8 


-HO. 


,17 


30. 


.92 


31. 


.69 


31. 


.74 


31. 


.76 


31. 


.94 


0. 


.60 




290 


172442 


.28-341128. 





-0.1 


21. 


.6 


-HO. 


.1 


-3, 


,3 


8.0 


+2 


,0 


-0.06 


64. 


,6 


-HO. 


,08 


30. 


.97 


31 


.55 


31 


.57 


31 


65 


31 


.76 


0. 


.61 


2MASS 


319 


172443 


.11-341215. 


9 


-0.2 


21 


.8 


-HO. 


.2 


-2 


,3 


5.6 


+4 


,4 


-0.1 


64. 


,6 


-HO. 


,1 


30. 


.87 


31 


.46 


31 


.49 


31 


.56 


31 


.71 


1 


.10 


USNO ()557-0528600-H2MASS 


331 


172443. 


.28-341243. 


.9 




21. 


.3 










1.4 








54. 


,5 






31. 


.08 


30. 


.76 


30. 


.78 


31. 


.25 


31. 


.43 


2. 


.61 


P; Pis24-2-HTYC- 7383-421-1 


335 


172443. 


.34-341137. 


.0 


-0.1 


22. 


.0 


-HO. 


,2 


-1, 


,4 


3.6 


-H2. 


,5 


-0.1 


64. 


,7 


-HO. 


,2 


30. 


.82 


31. 


.49 


31. 


.54 


31. 


.57 


31. 


.83 


1. 


.02 




344 


172443. 


.49-341156. 


.9 


-0.04 


22. 


.1 


-HO. 


.03 


-0.07 


1.2 


-HO. 


,07 


-0.05 


66. 


,2 


-HO. 


,05 


32. 


.19 


32. 


.26 


32. 


,36 


32. 


.53 


33. 


,11 


1. 


.57 


2T; Pis24-1-HTYC- 7383-7-1 


372 


172443 


.88-341139. 


.3 


-0.1 


21 


.8 


+0. 


.1 


-1 


,6 


5.2 


-H3 


,9 


-0.08 


54. 


,8 


+0. 


,09 


31. 


.07 


31 


.65 


31 


,68 


31 


.75 


31 


,91 


1 


.40 


2, "7 from Pis24-19 


379 


172443 


.95-341145. 


.6 


-0.03 


22 


.0 


+0. 


.06 




> 


10. 







-0.02 


55. 


,5 


+0. 


,05 


31. 


.55 


32 


.43 


32. 


,47 


32. 


.48 


32. 


,65 





.81 


Si 2MASS 


399 


172444. 


.37-341039. 


.8 


-0.09 


22 


.0 


+0. 


.09 


-1 


,1 


4.1 


-H2 


,0 


-0.09 


54. 


,9 


+0. 


,11 


31. 


.04 


31 


.74 


31. 


,79 


31. 


.82 


32. 


,05 


0. 


.95 


USNO 0558-0532335-H2MASS 


420 


172444. 


.72-341202. 


.6 


-0.03 


22 


.1 


-HO. 


,03 


-0.06 


0.7 


-HO. 


,06 


-0.08 


55. 


,9 


-HO. 


,13 


31. 


.73 


31. 


.42 


31. 


,55 


31. 


.90 


32. 


,83 


1. 


.16 


VT; Pis24-17 


464 


172445. 


.78-340939. 


.8 


-0.4 


21 


.1 


-HO. 


.3 


-0, 


,3 


2.0 


-HO. 


,5 


-0.08 


54. 


,4 


-HO. 


,08 


31. 


.03 


30 


,91 


30. 


,92 


31. 


.27 


31. 


,36 


1. 


.16 


Pis24-13 


471 


172446. 


.01-341407. 


.3 


-0.07 


21 


.9 


-HO. 


,07 


-1, 


,3 


5.7 


-H2. 


,4 


-0.06 


55. 


,2 


-HO. 


,06 


31. 


.35 


32 


,06 


32. 


,10 


32. 


.13 


32. 


,32 


0. 


.73 


2MASS 


472 


172446. 


.02-340942. 


.3 


-0.1 


22 


.0 


-HO. 


,1 


-0, 


,7 


2.6 


-HI. 


,2 


-0.1 


54. 


,7 


-HO. 


,2 


30. 


.78 


31 


,27 


31. 


,33 


31. 


.39 


31. 


,69 


0. 


.36 




493 


172446. 


.68-341233. 


.4 


-0.1 


22 


.0 


-HO. 


,1 


-1 


,2 


3.5 


-H2. 


,6 


-0.1 


54. 


,7 


+0. 


,2 


30. 


.88 


31 


,46 


31. 


,51 


31 


.56 


31 


,80 


0. 


.70 


USNO 0557-0528655-H2MASS 


519 


172447. 


.56-341048. 


.6 


-0.07 


22 


.1 


+0. 


,06 


-2 


,2 


7.5 


-H2 


,5 


-0.05 


55 


,2 


+0. 


,06 


31. 


.27 


32 


,13 


32. 


,18 


32. 


.19 


32. 


,38 


0. 


.93 


2MASS 


530 


172447. 


.95-341735. 


.3 


-0.4 


21 


.4 


+0. 


,3 


-0, 


,9 


3.1 


-HI. 


,6 


-0.1 


54. 


,4 


+0. 


,1 


30. 


.98 


31 


,17 


31 


,18 


31 


.38 


31 


,49 


0. 


.33 


2MASS 


577 


172450. 


.15-341243. 


.4 


-0.07 


22. 


.0 


-HO. 


,08 


-1, 


,0 


4.1 


-HI. 


,2 


-0.07 


55. 


,1 


-HO. 


,09 


31. 


.26 


31. 


,94 


31. 


,99 


32. 


.02 


32. 


,25 


1. 


.29 


USNO 0557-0528716-H2MASS 


612 


172453. 


.10-341525. 


.6 


-0.1 


22. 


.2 


-HO. 


,1 


-0, 


,4 


1.8 


-HO. 


,7 


-0.2 


64. 


,9 


-HO. 


,2 


30. 


.81 


31. 


,32 


31. 


,42 


31. 


.44 


31. 


,90 


0. 


.83 


USNO 0557-0528750-H2MASS 


634 


172454. 


.31-341209. 


.4 


-0.06 


22. 


.3 


-HO. 


,08 




> 


10. 







-0.05 


64. 


,8 


-HO. 


,09 


30. 


.53 


31 


,72 


31. 


,80 


31. 


.75 


31. 


,98 


1. 


.65 


S; 2MASS 


649 


172455. 


.55-341631. 


.3 


-0.10 


22 


.1 


-HO. 


,10 


-1 


,3 4,3 


-H2 


,5 


-0.1 


55. 


,1 


+0. 


,1 


31. 


.08 


31 


,86 


31 


,92 


31 


.93 


32. 


,18 


0. 


.59 


USNO 0557-0528794-H2MASS 


672 


172457. 


.87-341203. 


.9 


-0.04 


22 


.3 


-HO. 


,04 




> 


10. 







-0.03 


55. 


,3 


+0. 


,04 


31 


.17 


32 


,27 


32 


,34 


32 


.30 


32. 


,52 


0. 


.98 


S; 2MASS 


692 


172459. 


.74-340958. 


.7 


-0.1 


22. 


.5 


-HO. 


,1 


-0, 


,6 


2.2 


-HI. 


,1 


-0.2 


55. 


,2 


-HO. 


,3 


30. 


.46 


31 


,57 


31. 


,77 


31. 


.60 


32. 


,17 


0. 


.49 




746 


172608. 


.76-341115. 


2 




22. 


.0 










1.9 








65. 


,2 






31. 


.31 


31. 


,66 


31. 


,72 


31. 


.82 


32 


,18 


2. 


.19 


P; MSX6C G363.2261-H00.8296 


747 


172608. 


.83-341112. 


.5 


-0.03 


21 


.8 


-HO. 


,03 


-0, 


,1 


2.3 


-HO. 


,2 


-0.03 


66. 


,2 


-HO. 


,02 


32 


.46 


32. 


.78 


32 


,82 


32 


.96 


33 


,21 


1. 


.35 


VT; WR 93-HTYC7383-162-1 


753 


172610. 


.93-340843. 


2 


-0.07 


22. 


.5 


-HO. 


,10 




> 


10. 







-0.06 


64. 


,9 


-HO. 


,14 


30. 


.31 


31. 


.78 


31. 


,90 


31. 


.80 


32. 


,08 


0. 


.98 


S 


770 


172616. 


.76-341211. 


.3 


-0.1 


22. 


.0 


-HO. 


,1 


-1, 


,1 


3.8 


-H2. 


,2 


-0.1 


64. 


,7 


-HO. 


,1 


30. 


.92 


31. 


.52 


31. 


,67 


31. 


.62 


31. 


,84 


0. 


.74 


USNO 0557-0529090-H2MASS 



*A11 fits were "wabs(apec)" in XSPEC and assumed 0.3 X solar abundances. 

'^X-ray luminosities: s = soft band (0.5—2 keV); h = hard band (2—8 keV); t = total band (0.5—8 keV). Absorption-corrected luminosities are subscripted with a c. Luminosities were calculated assuming a 
distance of 2.56 kpc. 

'^2T means a two-tcmpcraturc fit was also performed; VT means a variable abundance fit was also performed; sec §6.1 for fit parameters. S means the spectral fit is only used as a spline fit for a rough 
estimate of luminosity; actual fit values are unreliable. P means the fit does not formally satisfy the convergence criterion, nevertheless the fit is acceptable to estimate the luminosity. The confidence intervals 
are missing since no error calculation is reported in XSPEC. Pis24 stars and other cataloged names obtained from VizieR are also shown. [/5iVO-The USNO-Bl.O Catalog (Monet et al. 2003); it TYC-The 



Tycho-2 Catalogue (Hog et al. 2000); MSX6C-MSX6C Infrared Point Source Catalog (Egan et al. 2003). 



Note. — Column 1: X-ray source number. Column 2: lAU designation. Columns 3-5: Estimated column density, plasma energy, and the plasma emission measure. Columns 6-10: Observed and corrected for 
absorption X-ray luminosities, obtained from our spectral analysis §3.4. Columns 11: Reduced fo^' the spectral fit. 



Table 4 

X-RAY Spectroscopy for Less Bright Sources: Thermal Plasma Fits 





Source 




Spectral Fit'' 






X- 


ray Luminosities'' 




Notes'^ 


Seq 


CXOU J 


log Nh 


kT 


log EM 


log Ls 


logLh 


log Lh,c 


log Lt 


log Lt,c 




# 




(cm-2) 


(keV) 


(cm^'^) 






(ergs s^-"-) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


3 


172359.11-341217.6 


22.0 


4.8 


53.9 


30.01 


30.73 


30.78 


30.81 


31.02 




8 


172405.42-341257.5 


-0.1 22.7 -HO.i 


>10 


54.7 


29.59 


31.60 


31.76 


31.61 


31.89 


S 


9 


172405.64-340710.0 


21.0 


-0.3 1.8 +0.3 


53.9 


30.62 


30.38 


30.39 


30.82 


30.86 




12 


172407.65-341751.9 


-0.2 22.3 +0.1 


-0.2 0.5 +0.1 


-0.2 55.0 +0.2 


30.36 


30.04 


30.25 


30.53 


31.93 




15 


172409.28-341224.3 


-1.5 22.5 +0.2 


-1.4 2.7 


-0.3 54.9 


30.30 


31.41 


31.58 


31.45 


31.92 




19 


172413.61-340921.8 


-0.5 21.5 +0.4 


-1.5 2.9 +14.8 


-0.3 54.0 +0.3 


30.53 


30.74 


30.76 


30.95 


31.08 




20 


172413.68-341451.8 


-0.4 22.2 +0.2 


-0.7 1.5 +2.6 


-0.6 54.4 +0.5 


30.24 


30.67 


30.79 


30.80 


31.36 




29 


172416.17-341526.6 


-0.2 22.4 +0.2 


-0.9 2.3 +4.3 


-0.4 54.7 +0.4 


30.30 


31.19 


31.34 


31.24 


31.73 




30 


172416.18-341311.3 


-0.9 21.6 +0.3 


-0.9 1.9 +3.5 


-0.3 54.1 +0.2 


30.49 


30.55 


30.58 


30.82 


31.05 




35 


172418.81-341204.7 


22.1 


1.3 


54.4 


30.43 


30.55 


30.65 


30.79 


31.37 





''AH fits were "wabs(apec)" in XSPEC and assumed 0.3 x solar abundances. 

''X-ray luminosities: s = soft band (0.5-2 keV); h = hard band (2-8 keV); t = total band (0.5-8 keV). Absorption-corrected luminosities are 

subscripted with a c. Luminosities were calculated assuming a distance of 2.56 kpc. 

'^S means the spectral fit is only used as a spline fit for a rough estimate of luminosity; actual fit values are unreliable. Some confidence intervals 
are missing since no error calculation is reported in XSPEC. 

Note. — Column 1: X-ray source number. Column 2; lAU designation. Columns 3-5; Estimated column density, plasma energy, and the plasma 
emission measure. Columns 6-10: Observed and corrected for absorption X-ray luminosities. The full table of 144 sources is provided in the electronic 
edition of the Journal. 



Table 5 

Stellar Counterparts and Classifications 





Source 




Counterparts & Class 




2MASS Coordinates 




2MASS Colors & Maf 


^nitudes 


Cautionary 


Flags 


Scq # 


CXOU J 


2MASS SIRIUS HST 


Flags 


RA{J2000) 


DEC(J2000) 


A0(") 


J - H 


H - Ki, 


Ks 


QFLAG 


Class 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


172353.98-340850.8 


X 






















2 


172358.42-340802.6 






R 


17 23 58.36 


-34 08 02.5 


0.87 


2.004 


1.261 


11.928 


AAA 


000 


3 


172359.11-341217.6 






R 


17 23 59.16 


-34 12 17.0 


0.80 


0.749 


0.487 


10.690 


AAA 


000 


4 


172401.20-340928.0 


%/ 




R 


17 24 01.32 


-34 09 27.3 


1.58 


0.744 


0.477 


12.926 


UAA 


000 


5 


172402.27-341110.4 


X 




H 


















(i 


172403.19-341402.2 






R 


17 24 03.22 


-34 14 02.1 


0.35 


0.990 


0.369 


11.313 


AAA 


000 


7 


172403.64-340634.7 


■J 




R 


17 24 03.56 


-34 06 36.2 


1.77 


2.092 


0.739 


13.505 


DBA 


ccc 


8 


172405.42-341257.5 


X 




H 


















9 


172405.64-340710.0 


v' 




F 


17 24 05.62 


-34 07 09.5 


0.60 


0.399 


0.172 


8.458 


AAA 


000 


10 


172406.70-341305.0 


V 




HR 


17 24 06.64 


-34 13 04.9 


0.74 


2.319 


1.208 


13.316 


UCA 


000 



Note. — Column 1: X-ray source number. Column 2: lAU designation. Columns 3-5: existence of ONIR counterparts; -y^-counterpart. X-no counterpart, ■ - - -source is out of the 
field of view of the available images. Column 6: X-ray source classifications. See §5.5 for descriptions of the class designations. Column 7-9: 2MASS (J2000) coordinates and the offset 
between X-ray position and near-lR position in arcsec. Column 10-12: 2MASS J — H,H — Kg, and Kg measurements. Column 13: From Explanatory Supplement to the 2MASS All Sky 
Data Release (Cutri et al. 2003). 2MASS JHKs photometry quality flags with three characters: A — very high significance detection (> 10 SNR); B — high significance detection (> 7 
SNR,); C — moderate significance detection (> 5 SNR); D — low significance detection; E — point spread fitting poor; F — reliable photometric errors not available; U — upper limit 
on magnitude (source not detected); X = source is detected but no valid brightness estimate is available. Column 14: From Explanatory Supplement to the 2MASS All Sky Data Release 
(Cutri ct al. 20t)3). 2MASS JHKs confusion and contamination flags with three characters: = source unaffected by artifacts; b = bandmerge confusion (possible multiple source); c 
— photometric confusion from nearby star; d — diffraction spike confusion from nearby star; p = persistence contamination from nearby star; s = electronic stripe from nearby star. 
The full table of 779 Chandra sources is available in the electronic edition of the Journal. 



Table 6 

Known O and Early B stars in Pismis 24 



Star ONIR photometry X-ray properties 



ID" 


SpecType'' 


«,T2000 


'5.T2000 ^ 




U - 




V - J"^ 


J - H'^ 


H - Ks'^ 


if.,'' 




Scq # 




logJV^s 


logibol'' 


Notes'' 


Pis 24-1 


03If* 


17 24 43.49 


-34 11 57.0 


10.43 


0.40 


1.45 


3.70 


0.55 


0.28 


5.89 


6.16 


#344 


33.11 


22.1 


39.60 


HDE 319718=[N78] 35 


Pis 24-17 


03III(f*) 


17 24 44.73 


-34 12 02.7 


11.84 




1.49 


4.02 


0.54 


0.31 


6.98 


6.57 


#420 


32.83 


22.1 


39.09 


[N78] 57 


Pis 24-2 


06.6V((f)) 


17 24 43.20 


-34 12 43.5 


11.95 


0.32 


1.41 


3.48 


0.48 


0.31 


7.69 


5.81 


#331 


31.43 


21.3 


38.81 


Pis 24-13 


06.5V((f)) 


17 24 45.78 


-34 09 39.9 


12.73 


0.11 


1.48 


3.92 


0.57 


0.29 


7.95 


6.45 


#464 


31.36 


21.1 


38.61 


[N78] 36 


Pis 24-16 


07.5V 


17 24 44.45 


-34 11 58.8 


13.02 




1.60 


4.57 


0.59 


0.49 


7.37 


7.35 


#407 


31.90 


22.0 


38.60 


[N781 58 


Pis 24-3 


08V 


17 24 42.30 


-34 13 21.5 


12.75 


0.24 


1.41 


3.46 


0.52 


0.30 


8.47 


5.79 


#293 


31.17 


21.0 


38.57 




Pis 24-15 


08V 


17 24 28.86 


-34 14 50.3 


12.32 


0.14 


1.27 


3.29 


0.45 


0.24 


8.35 


5.54 


#77 


32.04 


21.5 


38.63 


[N78] 46 


Pis 24-10 


09V 


17 24 3fi.04 


-34 14 00.5 


13.02 


0.40 


1.40 


3.51 


0.49 


0.28 


8.75 


5.84 


#151 


31.96 


22.1 


38.36 


WR 93 


WC7+07-9 


17 25 08.85 


34 11 12.5 


10.68 


0.70 


1.42 


3.65 


0.50 


0.67 


5.87 


5.86 


#747 


33.21 


21.8 


>39.97 


HD 167604= [N78] 62 


Pis 24-18 


B0.5V 


17 24 43.29 


-34 11 41.9 


13.97 




1.48 


3.97 


0.68 


0.34 


8.97 


6.42 


#332 


31.56 


22.1 


38.26 




Pis 24-12 


BIV 


17 24 42.27 


-34 11 41.2 


13.88 


0.38 


1.47 


3.38 


0.51 


0.21 


9.78 


5.56 


No 






38.05 




Pis 24-19 


BIV 


17 24 43.69 


-34 11 40.7 


14.43 




1.39 


3.76 


0.59 


0.47 


9.61 


6.08 


No 






38.04 




Pis 24-8' 


09? 


17 24 38.81 


-34 14 58.2 


12.98 


0.48 


1.44 


3.67 


0.59 


0.21 


8.52 


5.90 


#194 


31.29 


21.9 






Pis 24-7' 


09? 


17 24 47.81 


-34 15 16.5 


13.46 


0.58 


1.68 


4.10 


0.61 


0.33 


8.42 


6.61 


#527 


30.28 


22.2 






Pis 24-4' 


O9-B0? 


17 24 40.39 


-34 12 05.9 


13.93 


0.53 


1.43 


3.88 


0.60 


0.32 


9.14 


6.49 


#225 


30.46 


22.0 






Pis 24-9' 


09-B1? 


17 24 39.29 


-34 15 26.4 


14.26 


0.40 


1.40 


3.48 


0.49 


0.29 


10.00 


5.81 


No 










Pis 24-11' 


09-B1? 


17 24 34.68 


-34 13 17.1 


14.53 


0.30:: 


1.67 


3.79 


0.66 


0.30 


9.79 


7.03 


#140 


32.20 


21.9 







^Identifications from Moffat & Vogt (1973) and Massey et al. (2001). 

'^Spectral types from Massey et al. (2001) and Walborn et al. (2002). Optical photometry from Massey et al. (2001). Bolometric luminosities from Massey et al. (2001) are recalculated using bolometric 
corrections from Heap et al. (2006) instead of Vacca et al. (1996). 

*^Units of RA are hours, minutes and seconds; units of DEC are degrees, arcminutes, and arcseconds. Coordinates are from the Catalogue of Stellar Spectral Classifications (Skiff 2005) and Massey et al. 
(2001). 

■^Infrared photometry from 2MASS. 

""Visual extinction from Bohigas et al. (2004) (Ay ^ 1.39E{V - J)). 
^ACIS source number for X-ray detections. 

^Lx is absorption corrected full-band luminosity ^ (ergs s""*^). is X-ray derived column density (cm~^). 

^Other identifications from Neckel (1978). 

^Not included in Massey et al. (2001) spectroscopic classification sample. Tentative spectral types are inferred from color-magnitude diagram by de-reddening 2MASS photometric colors to the 1 Myr 
pre-MS isochrone. 



Table 7 

x-ray and ir luminous stars as candidate o stars 



Seq# 


CXOU J 


J-H 


H-K. 




logLx 


log Nh 


kT 


Notes 








(a) 


Stars with log Lx > 32.0 






140 


172434.79-341318.0 


0.660 


0.296 


9.789 


32.04 


21.9 


7.0 


Pis 24-11 


158 


172436.63-341550.7 


1.194 


0.504 


12.270 


32.01 


22.0 


2.9 


V 


227 


172440.67-341403.8 


1.203 


0.176 


14.324 


32.06 


22.1 


7.6 


flaring 


379 


172443.95-341145.6 


-0.572: 


0.651: 


12.298 


32.65 


22.0 


10.0: 


HST 


399 


172444.37-341039.8 


1.063 


0.957 


10.972 


32.05 


22.0 


4.1 


HST 


471 


172446.01-341407.3 


1.636 


0.802 


10.674 


32.32 


21.9 


5.7 




519 


172447.56-341048.6 


1.466 


0.706 


10.167 


32.38 


22.1 


7.5 


V;HST 


577 


172450.15-341243.4 


1.240 


0.509 


10.788 


32.25 


22.0 


4.1 




649 


172455.55-341631.3 


1.286 


0.510 


12.064 


32.18 


22.1 


4.3 


V 


672 


172457.87-341203.9 


1.331 


0.576 


10.256 


32.52 


22.3 


10.0: 


flaring 


692 


172459.74-340958.7 








32.17 


22.5 


2.2 


faint in SIRIUS 


746 


172508.76-341115.2 








32.18 


22.0: 


1.9: 


close to WR star 


753 


172510.93-340843.2 








32.08 


22.0 


3.8 


V; IR nebulosity 


(b) Stars with K < 10 


9 


172405.64-340710.0 


0.399 


0.172 


8.458 


31.20 


21.5 


0.7 


IRAS 17207-3404 


19 


172413.61-341656.5 


2.006 


0.991 


5.334 


31.0 


21.9 


2.8 


MSX6C G353.0401-I-00.9326 


120 


172433.46-341344.9 


0.602 


0.333 


9.649 


31.29 


21.7 


10.0: 




194 


172438.90-341459.0 


0.587 


0.208 


8.519 


31.29 


21.9 


1.7 


Pis24-8 


225 


172440.50-341206.3 


0.595 


0.315 


9.138 








Pis24-4 


312 


172442.88-340911.6 


0.931 


0.502 


9.231 










527 


172447.87-341517.0 


0.608 


0.327 


8.424 








Pis24^7 


539 


172448.44^341743.7 


0.468 


0.270 


8.026 


31.14 


21.5 


0.9 




643 


172455.09-341111.5 


0.798 


0.427 


8.866 


32.15 


22.4 


0.6 




654 


172455.85-341234.0 


1.297 


0.952 


8.389 










694 


172459.86-341610.8 


1.645 


1.304 


8.137 


31.34 


22.1 


2.5 


MSX6C G353. 1413-1-00.8083 



Note. — Column 1: X-ray source number. Column 2: lAU designation. Columns 3-5: 2MASS colors and Kg magnitudes. 
Columns 6-8: X-ray luminosities (0.5-8.0 keV) corrected for absorption, derived column densities, and thermal plasma 
energy obtained from our spectral analysis. Columns 9: Related notes: ACIS source number; V=Variable lightcurves; 
HST=HST counterparts found; Other catalog names. Stars in the top panel are selected from X-ray criteria, and those 
in the bottom panel are selected from infrared criteria. Missing logLx, logATjj, and kT values mean that the source has 
< 20 net counts. 



Table 8 

X-RAY Selected Candidate Intermediate- mass Stars 



Seq# 


CXOU J 


J 


H 


Ks 












(mag) 


(mag) 


(mag) 


(mag) 


(Mq) 


(mag) 


2 


172358.42-340802.6 


14.99 


13.16 


11.93 


18 


13 




3 


172359.11-341217.6 


11.86 


11.16 


10.69 


8 


14 




6 


172403.19-341402.2 


12.60 


11.67 


11.31 


10 


13 




12 


172407.65-341751.9 


14.32 


13.03 


12.46 


8 


3 




18 


172413.61-340921.8 


13.51 


12.65 


12.34 


4 


3 




20 


172413.68-341451.8 


13.12 


12.13 


11.77 


7 


4 




29 


172416.17-341526.6 


13.88 


12.28 


11.45 


16 


16 




35 


172418.81-341204.7 


14.80 


13.64 


13.04 


6 


2 




38 


172420.25-341310.7 


13.90 


12.54 


11.88 


14 


12 




45 


172421.75-341539.4 


15.40 


13.55 


12.82 


15 


4 




156 


172436.45-341236.7 


13.30 


12.32 


11.57 


5 


3 


V 



Note. — Column 1: X-ray source number. Columns 2: lAU designation. Columns 3-5: 
2MASS JHKs-hand magnitudes. Columns 6-7: Visual extinction and mass estimated 
from de-reddened location along the standard interstellar reddening vector to the 1 
Myr prc-MS isochronc in the 2MASS color-magnitude diagram, assuming a distance 
d ~ 2.56 kpc. Column 8: Flag of significant if -band excess. The full table of 100 
Chandra sources is available in the electronic edition of the Journal. 

^These quantities assume a cluster age of 1 Myr and d ~ 2.56 kpc. For early B stars, 
Av and M (estimated visual extinction and mass) are especially dependent on the age 
and distance assumptions (Figure 7). 

''x/ marks significant Jf-band excess derived from high quality 2MASS photometry 
where near-IR color excess E{H — K) is larger than cr{H — K), the uncertainty in H — K 
color index. 



Table 9 

X-RAY Selected Heavily Obscured Sources 



Seq # CXOU J Criteria'' logN„ {cm-'^)'° (E) (keV)" Notes'^ 



8 


172405.42-341257.5 


I;II 


22.7 


3.9 




15 


172409.28-341224.3 


I 


22.5 


3.1 


V;2MASS 


61 


172426.47-341850.9 


I 


22.5 


3.3 




142 


172134.95-340533.3 


I;II 


23.2 


4.7 




569 


172449.82-341455.6 


I 


22.4 


3.4 




598 


172451.90-341331.3 


I;II 


22.6 


3.1 


2MASS 


599 


172451.95-341513.0 


II 


22.6 


2.0 


2MASS 


634 


172454.31-341209.4 


I 


22.3 


3.1 


VV;2MASS 


664 


172456.97-341651.8 


II 


22.6 


2.1 


2MASS 


692 


172459.74-340958.7 


II 


22.5 


2.9 




693 


172459.85-341302.9 


II 


22.5 


2.8 


2MASS 


702 


172501.09-341404.6 


I;II 


22.7 


3.0 




713 


172502.54-341549.7 


I;II 


22.7 


4.1 




728 


172504.93-340410.5 


I;II 


22.6 


3.9 


e 


753 


172510.93-340843.2 


I;II 


22.5 


3.5 


VV;f 


775 


172519.23-341356.6 


I;II 


22.8 


3.7 





^Selection criteria: 1-{E) > 3.0 keM; 11.-Nh,x > 22.5 cm ^. The values shown here 
are rounded. 

''Column density derived from spectral fitting. 

''Median photon energy of the X-ray source. 

■^Related notes: V=variable Ughtcurves; VV=flaring; 2MASS=2MASS counterparts 
found. 

"^1^8X60 G353.3175-f00.9075 
''MSX6C G353.2651-1-00.8443 
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Fig. 1. — Multiwavelength images of the NGC 6357 region from long to short wavelengths, (a) 
1.5° X 1.5° view of 6 cm radio continuum contours from Parkes survey (Haynes et al. 1978) with 
the three Hii regions labeled, (b) The same region shown as a 3-color composite mid-infrared 
image from the MSX survey: red is 8.28 /xm, green is 12.13 //m, and blue is 21.3 /im. (c) 20' x 20' 
K-band image from the 2MASS survey showing the Hii region G353.2+0.9, dark clouds, and the 
stellar field. The box shows the SIRIUS FOV. (d) 6' x 6' JHK composite image from the 
SIRIUS observation. J is shown in blue, H is green, and K is red. (e) 20' x 20' R-band image 
from the Digitized Sky Survey emphasizing Ha emission. The box shows HST/WFPC2 coverage. 
The brightest members WR 93, Pis 24-1, and Pis 24-17 are marked, (f) An unpublished X-ray 
image from a short ROSAT PSPC observation. The three brightest regions are associated with 
the NGC 6357 Hii regions. The ACIS-I FOV is shown as a box, with the off-axis S2+S3 chips 
represented as a rectangle to the east. Note that S3 captures G353.2+0.7 at its southern edge. 
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Fig. 2.— (a) The 17' x IT ACIS-I image (0.5-7.0 keV) binned by 2" x 2". The grey polygons 
show the source extraction regions that are matched to the local Chandra point spread 
function. The black contour in the center outlines the ionization front seen in the optical 
image, (b) The ACIS full-field adaptively smoothed image showing all four I-chips and two 
S-chips. Red represents the soft band (0.5-2.0 keV) X-ray emission and blue represents the 
hard band (2.0-7.0 keV) X-ray emission. The cyan and white polygons outline the extraction 
regions for possible diffuse emission in the central region, and in a larger area, respectively. 
The yellow polygon outhnes the region used for diffuse background estimation. See §7. (c) 
Zooming-in on the inner ~ 4' x 4' of the ACIS-I data, binned by 0.25" x 0.25", then adaptively 
smoothed. 
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Fig. 2. — Continued. 
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Fig. 3.— Examples of close X-ray source pairs (ACIS #313, 317; #402, 406; #427, 428; shown as 
dashed circles) recovered from image reconstruction. Each panel includes: 6" x 6" raw ACIS image 
binned by 0.25" (top left), reconstructed image (top right), HST image (bottom left), and SIRIUS 
image (bottom right). The white stripe in the HST image is the diffraction spike from a nearby 
bright star. 




Fig. 4. — (a) HST/WFPC2 image (F814W) revealing details of the interface between the 
massive stars, the H ll region, and the molecular cloud; (b) zoom in on the elephant trunk 
with 2MASS near-IR sources (red circles) and X-ray sources (cyan polygons) overlaid. The 
cyan polygons are not circles; they represent the X-ray source extraction polygons based on 
the PSF. 
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Fig. 5. — A cartoon illustration of the source classification in §4.5: i?-Reliable X-ray sources 
with optical and NIR stellar identifications; C-Clustered X-ray sources with off-axis angle 
< 5.0'; i/-Hard X-ray sources with {E) > 3.0 keV. The area of each circle is proportional to 
the number of sources in the class. Group U (unreliable X-ray detections) and F (foreground 
field stars) are not shown here. 
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Fig. 6. — Color-color diagram of X-ray stars with high quality 2MASS counterparts. Open 
circles and filled circles represent sources without and with significant X-band excess {E{H — 
K) > 1(7 uncertainty in H — K color), respectively. Crosses mark the known members of 
the Pismis 24 cluster from Massey et al. (2001). ACIS source numbers and error bars are 
given for the counterparts with significant 7^-band excess. The interesting source ^^19, the 
Pis24 member #331, and the foreground star #773 are also labeled to guide the reader. The 
solid and long-dash lines denote the locus of main sequence stars and giants from Bessell & 
Brett (1988), respectively. The dash dotted line is the locus for T Tauri stars from Meyer 
et al. (1997) and the dash triple-dotted line is the locus for HAeBe stars from Lada & 
Adams (1992). The dash lines represent the standard reddening vector for Ay = 20 mag, 
with asterisks marking every Ay = 5 mag. Many additional X-ray stars with lower quality 
2MASS photometry are omitted. 
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Fig. 7. — Infrared J—H vs. J color- magnitude diagram. The solid line is the 1 Myr isochrone 
for pre-MS stars from Baraffc ct al. (1998) and Siess et al. (2000). A 0.3 Myr isochrone (Siess 
et al. 2000) is also shown (dotted line) to demonstrate how the estimated masses differ if we 
adopt a younger age for the cluster. The dash dotted line marks the location of Zero Age 
Main Sequence (ZAMS) stars. Dashed lines represent the standard reddening vector with 
asterisks marking every Ay — 5 mag; the corresponding stellar masses are marked. ACIS 
numbers for sources with significant i^-excess are shown. The interesting source #19, the 
Pis24 member #331, and the foreground star #773 are labeled to guide the reader. Error 
bars are shown for all sources with color uncertainties > 0.1 mag. This diagram uses the 
same source sample and symbols as Figure 6. 
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Fig. 8. — Infrared H — K vs. K color-magnitude diagram. The solid line is the 1 Myr 
isochrone from Baraffe et al. (1998) and Siess et al. (2000). The dash dotted line marks the 
location of ZAMS stars. Dashed lines represent the standard reddening vector with asterisks 

marking every Ay = 5 mag; the corresponding stellar masses are marked. ACIS numbers for 
sources with significant i^T-cxccss arc shown. The interesting source 7^19, the Pis24 member 
#331, and the foreground star #773 are labeled to guide the reader. Error bars are shown 
for all sources with color uncertainties > 0.1 mag. This diagram uses the same source sample 
and symbols as Figure 6. 
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Fig. 9. — (a): X-ray luminosity function (XLF) constructed from the observed (not cor- 
rected for absorption) hard band (2.0-8.0 keV) X-ray luminosity for the lightly obscured 
population in NGC 6357 (solid line), Orion (dashed line, Feigelson et al. 2005), and Cep B 
(dotted line, Getman et al. 2006). (b): XLF using for the heavily obscured population 
in NGC 6357 and Cep B compared to the XLF for the lightly obscured population in Orion, 
(c): XLF using absorption corrected hard band (2.0-8.0 keV) X-ray luminosity L^ f. for the 
lightly obscured population in NGC 6357, Orion, and Cep B. (d): XLF using L^ c for the 
heavily obscured population in NGC 6357 and Cep B compared to the XLF for the lightly 
obscured population in Orion. 
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Fig. 10. — The spatial distribution of the X-ray stellar sources. Members with ONIR counter- 
parts are plotted as diamonds, and likely members without ONIR identifications are shown 
with crosses. Obscured sources {{E) > 3.0 keV) are plotted in blue and less obscured sources 
{{E) < 3.0 keV) are plotted in red. The large scale optical nebulosity contour from the DSS 
image traces the mid-IR ring-like morphology; it is shown in black together with the outlines 
of the ionization front in yellow and the CO "South-Eastern Complex" in green. 
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Fig. 11. — Projected stellar surface density (number of stars arcmin"^) maps for (a) the 
unobscured population and (b) the obscured population, shown in logarithmic scale using a 
0'.5 radius smoothing kernel. The highest concentration of stars coincides with the core of the 
cluster, where the massive O stars are located. A density enhancement, or a subcluster is seen 
centered at RA=17.414 (hours), Dec=-34.24 (degrees) in (b). There is a small overdensity 
around the WR star, at RA=17.419 (hours), Dec=-34.19 (degrees) in (a). 




Fig. 12. — Radial density profiles for X-ray stars in the Orion Nebula Cluster (COUP 
sample, Getman et al. 2005; thin line) and for Pismis 24 (thick line). COUP sample is 
spatially complete in the central 1 pc of Orion Nebula Cluster. 
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Fig. 13. — The distribution of X-ray sources overlaid on the central 50" SIRIUS K image. 
The ONIR identified sources (circles), and unidentified sources (crosses) are shown, with the 
known massive members highlighted in white. 
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Fig. 14. — X-ray spectra of the 03 stars and the brightest source WR93. The left three 
panels give apec model fits for Pis24-1 {Nh = 1.1 x 10^^ cm-^ and kT = 1.2 keV), Pis24-17 
{Nh = 1.2 X 10^2 cm-2 and kT = 0.7 keV), and WR 93 {Nh = 6.9 x lO^i cm'^ and kT = 2.3 
keV). Top right shows an improved fit with a two temperature thermal plasma fit for Pis24-1 
{Nh = 1.3 X 10^2 cm-2, kTi = 0.5 keV and kT2 = 1.7 keV). Middle right gives a slightly 
improved fit for Pis24-17 with a variable abundance thermal plasma fit {Nh = 1.4 x 10^^ 
cm"2 and kT — 0.7 keV). Bottom right presents an improved fit for WR 93 with the vapec 
model {Nh = 7.4 x lO^^ cm-^ and kT = 2.4 keV). 
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Fig. 15. — The Lx vs. L^oi relation for the X-ray detected OB stars in Pismis 24 (diamonds) 
compared to COUP massive stars with strong winds (crosses) and weak winds (squares). 
Four COUP stars with weak winds have upper hmits in Lx reported in Stelzer et al. (2005). 
The L\)oi for WR 93 is a lower limit due to uncertain bolometric correction for late-type WC 
stars. The dashed line represents the canonical relation Xo'giLxli-'hoi) ^ 




Fig. 16. — Photon arrival diagram, lightcurve (2 ks bins) and X-ray spectrum of source #672 
(CXOU J172457.87-341203.9), with a powerful flare detected in this ACTS observation. The 
best spectral fit is obtained with Nh = 1.7 x 10^^ cm~^ and kT > 10 keV. 
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Fig. 17. — The spatial distribution of deeply embedded sources. Sources that have {E) > 3.0 
keV are plotted with crosses and sources that have log Nh,x > 22.5 are plotted with circles. 
Diamonds represent sources that satisfy both criteria. The background DSS image is inverted 
so that brighter regions appear darker. 



